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Abstract. We study many-body properties of quantum harmonic oscillator lattices with 
disorder. A sufficient condition for dynamical localization, expressed as a zero-velocity Lieb- 
Robinson bound, is formulated in terms of the decay of the eigenfunction correlators for 
an effective one-particle Hamiltonian. We show how state-of-the-art techniques for proving 
I Anderson localization can be used to prove that these properties hold in a number of standard 

■ models. We also derive bounds on the static and dynamic correlation functions at both zero 

and positive temperature in terms of one-particle eigenfunction correlators. In particular, 
we show that static correlations decay exponentially fast if the corresponding effective one- 
particle Hamiltonian exhibits localization at low energies, regardless of whether there is a 
gap in the spectrum above the ground state or not. Our results apply to finite as well as to 
CO ' infinite oscillator systems. The eigenfunction correlators that appear are more general than 

those previously studied in the literature. In particular, we must allow for functions of the 
Hamiltonian that have a singularity at the bottom of the spectrum. We prove exponential 
. bounds for such correlators for some of the standard models. 

1. Introduction 

Oscillator lattice systems are the standard model for the vibrational degrees of freedom, 
known as phonons, in crystal lattices. These phonons interact with the other degrees of 
freedom, such as spins and electrons, in ways that often significantly modify their behavior. 
The two-point correlation functions in the ground state and thermal states of the lattice, e.g. 
CN ' (ixQy) an d (PxPy), and their time-dependent analogues, are key quantities to understand the 

qq ■ role of the lattice variables in the structure and transport properties of materials. 

For a pure regular crystal in the harmonic approximation, finding the correlation functions 
is a standard exercise. It has also long been understood, however, that to describe a basic 
phenomenon such as heat conduction, it is important to go beyond the pure harmonic lattice 
and to consider anharmonicities [2H] and randomness [3B]. Anharmonicities lead to non-linear 
phonon interactions, as studied, e.g., in |40j . In this work we are interested in the effects of 
randomness, which could e.g. occur in a crystalline material layered on a substrate. 

The main effect of randomness is that it may lead to localization of the phonons at certain 
energies, which means that a portion of the energy spectrum consists of eigenvectors that 
decay exponentially fast at large distances. When localization occurs at all energies the 
physical effect is dynamical localization, i.e., the absence of propagation of disturbances in 
the system. This phenomenon has been studied extensively in one-particle models such as 
the Anderson model, see [HI [23 ESI H2] for surveys and more references. The mathematical 
study of localization in systems with many degrees of freedom has only recently begun, e.g. 

pin si si]. 
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There is no consensus on how to mathematically quantify dynamical localization in the 
many-body context. With this as motivation, we will presently, as indicated above, consider 
random oscillator systems. For such systems, much can be analyzed directly and in fact, we 
will prove two main estimates which, under certain assumptions, suggest that the system is 
dynamically localized. First, we give conditions under which the models satisfy a zero velocity 
Lieb- Robinson bound. Next, we prove that time-dependent correlations (both ground state 
and thermal) decay exponentially in space with no gap assumption. Below, we briefly describe 
both of these results. 

Our first results demonstrate a condition under which the random oscillator systems we 
consider satisfy a strong form of dynamical localization; namely, a zero velocity Lieb- Robinson 
bound (in disorder average). Generally, Lieb- Robinson bounds are deterministic and show 
that the Heisenberg evolution (corresponding to a Hamiltonian comprised of sufficiently short 
range interactions) of a local observable remains essentially local. More precisely, they prove 
that given an observable initially supported in a finite set X C the dynamic evolution of 
this observable to time t has, up to exponentially small corrections, support contained in a 
ball centered at X of radius v\t\. The number v > is called the group velocity or velocity of 
propagation for the system. Results of this type were originally proven by Lieb and Robinson 
[28] in the context of quantum spin systems, and they have recently been generalized (351 123] , 
in particular to the setting of oscillators systems, see e.g. [Ml EES1 E|- 

Our results here, see e.g. Theorems 13.31 and 13 . 8 1 for finite oscillator systems and Theorem 17. II 
for an extension to infinite systems, give assumptions under which random oscillator systems 
satisfy a Lieb-Robinson bound, on average, uniformly in time, i.e. a bound with velocity of 
propagation. Similar results were recently obtained in [21] for the random XY chain. 

The condition for the above results to hold constitutes our main new technical input. It 
is described in terms of estimates on certain singular eigenfunction correlators corresponding 
to the effective one-particle Hamiltonian. In the single particle theory of random systems, 
eigenfunction correlators are a well-known tool for deriving results on dynamical localization, 
see e.g. [H HI 02]. Here the new insight is that, for the effective models corresponding to 
random oscillator systems, singular eigenfunction correlators need to be estimated, and we 
demonstrate, see specifically Appendices \A\ and El that the desired bounds can indeed be 
proven in certain well-studied systems. 

Our next results concern correlation decay. It is well-known that, deterministically, gapped 
ground states satisfy exponential clustering, i.e. the correlations of observables in gapped 
ground states decay exponentially in the distance between the supports of the observables, 
see [351 [23] (spin systems) and [391 [151 [33] (oscillator systems) for proofs of this fact. 

For random oscillator systems, we give conditions under which ground state correlations 
satisfy exponential clustering, on average; even in cases where the gap closes. Again, the 
conditions for these results also involve estimating certain singular eigenfunction correlators. 
We note that such a situation, one for which exponential clustering holds without the existence 
of a uniform gap, can be referred to as a regime with a mobility gap, and this has been discussed 
in the literature, see e.g. [22]. Specifically, our results show that, under certain assumptions, 
exponential clustering holds, on average, for both static and dynamic correlations in ground 
states and also in thermal states. 

We now give a detailed outline of the contents of this paper. In Section [21 we introduce a 
class of oscillator systems depending explicitly on three sets of parameters: the masses, the 
spring constants, and the interaction strengths. To establish our notation, we briefly review 
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the well-known procedure by which these models can be diagonalized and comment on some 
important consequences. Here, as well as in the following Sections E] to [HI we consider finite 
(but arbitrarily large) oscillator systems. 

In Section [31 we prove our strongest form of localization for the dynamics corresponding 
to random oscillator systems. This localization is expressed in terms of a zero-velocity Lieb- 
Robinson bound, in average, and it follows from an appropriate exponential decay of the 
eigenfunction correlators corresponding to the effective one-particle Hamiltonian. We show 
that these results apply to two classes of observables: in Section 13.11 we consider Weyl opera- 
tors, while in Section 13.21 our bounds for the positions and momenta are discussed. 

In Sections @] and we present our results on dynamic and static correlation decay, again 
for both Weyl operators as well as the positions and momenta. Section H] discusses correlations 
in the ground state, whereas Section demonstrates our results for thermal states. 

In Section we focus on our main application — the case in which the effective one- 
particle Hamiltonian becomes the Anderson model. Here we illustrate how the bounds on the 
many-body system which we derived in Sections [31 HI and [5] depend on certain well-studied 
localization regimes for the Anderson model, e.g. large disorder and one- dimension. We also 
describe what our results show in cases where dynamical localization has only been proven at 
the bottom of the spectrum. 

The extension of our results from finite to infinite oscillator systems is the content of Sec- 
tion [TJ We focus on the case where the effective one-particle Hamiltonian is the infinite volume 
Anderson model. Again we obtain zero-velocity Lieb-Robinson bounds as well as exponential 
decay of ground state and thermal state correlations. 

In Section [8] we provide a brief non-technical summary of our new results. 

Two appendices contain results on Anderson localization for the effective one-particle op- 
erators associated with the oscillator system Hamiltonians, separate for finite volume (Ap- 
pendix [AJ and infinite volume (Appendix [B]) . Many of these results, in particular those on 
singular eigenfunction correlators, have not previously appeared in the literature. Thus we 
will present them with detailed proofs. We think that these results are of independent interest 
in the theory of Anderson localization. 

A number of the conditional results in this paper, see specifically those in Sections [3J HJ 
and apply to much more general systems (oscillator systems over more general graphs, 
more general interactions, see Remark 12.21 below). However, for the sake of clarity, we have 
compromised obvious extensions to illustrate our main techniques. Moreover, our primary 
application to the case of Anderson-type single particle Hamiltonians is of the form discussed, 
and more general results on dynamical localization for the corresponding single particle models 
are generally not known. 

A related topic of great current interest are the so-called Area Law bounds for the entropy of 
entanglement in the ground state of multi-dimensional systems. The techniques of this paper 
can be used to obtain such a bound for the random oscillator lattices considered here, under 
the assumption that there is a mobility gap (but not necessarily a gap) above the ground 
state. We will report on that result in a separate article [36] . 
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2. Systems of coupled harmonic quantum oscillators 

We will consider oscillator models defined over Z d , d > 1 an integer. These we describe in 
terms of three sequences of parameters: the masses {m x }, the spring constants {k x }, and the 
couplings {X X)V }. Throughout this work, we will assume the following. 

Assumption A.l: For each x £ Z d , m x > 0, k x > 0, and for each two element set 
{x,y} C Z d of nearest neighbor sites (e.g. \x — y\ = 1, where | ■ | will denote the 1-norm), 
^x, y > 0. Moreover, we assume there are numbers m min , m max , /c max , and A max such that 
< 7Tt m i n < vn x < ?n max , < k x < /c max , and < \ x . y < A max . 

In terms of these sequences, finite-volume harmonic Hamiltonians are defined as follows. 
Fix an integer L > 1 and a cubic set = [-L, L] d D Z d . Set 

(2.1) ^=E(^ + y^)+ E w&-?«) a 

X£A L V X / {x,y}CA L 

to be the harmonic Hamiltonian on A^ which acts on the Hilbert space 

(2.2) H L = L 2 (R A \dq) = (g) L 2 (R,dq x ), 

x&K L 

with q = {q x ) x( zA L denoting the space variables. Note here that we use subsets {x,y} C 
rather than pairs (x,y) G A^x A^, meaning that {{x,y} : \x — y\ — 1} labels the undirected 
edges of the next-neighbor graph associated with A^. Moreover, we also use q x to denote 
the position operators, meaning the multiplication operators by q x , and by p x = —id/dq x we 
denote the momentum operators. In this case, Hl describes a <i-dimensional system of one- 
dimensional harmonic oscillators that are coupled by nearest neighbor quadratic interactions. 
We consider finite systems of side- length 2L + 1, but will be interested in results (and bounds) 
which hold uniformly in the size of the system. This is one reason why we assume the uniform 
bounds on the parameters found in Assumption A.l. 

By standard results, e.g. [37], each q x , p x , and also Hl is essentially self-adjoint on C7^°(lR Ai ), 
the smooth functions of compact support, as well as on the Schwarz space of rapidly decreasing 
functions over R Al . From now on we will use to denote the self-adjoint closure of this 
operator. Clearly, is non-negative. 

As is well-known, it is possible to describe the eigenvalues and eigenfunctions of quite 
explicitly. This is due to the fact that the harmonic Hamiltonian Hl can be reduced to a 
system of uncoupled harmonic oscillators (or free Boson system) via a quantum canonical (or 
Bogoliubov) transformation, see, e.g., [18] or [10] for background. We review this procedure 
in order to introduce some relevant notation. To start, we express Hl in the form 

(2.3) HL = (q T ,p T )M L ( q p y 
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where 

(2.4) ^=(V I 

Here we view q — (q x ) and p = (p x ) as column vectors indexed by x G A^, with the transposes 
q T and p T denoting the corresponding row vectors. With this convention (12. 3 P can be under- 
stood in the sense of standard matrix multiplication. By /i we denote the diagonal matrix 
fi XtV = 2^~$x,y, which we will generally think of as a multiplication operator on £ 2 (Al) (omit- 
ting the L-dependence in our notation). Thus we may write f/, X)y = (5 x ,fiS y ), with {S x } x£ \ L 
the canonical basis in £ 2 (Al). 

The matrix-elements of the operator /io,L in £ 2 {Al) are given by 

(2-5) (5 X , h , L S y ) = <j -\ x , y , if \x-y\ = 1, 

0, else. 



The associated form is given by 

(2.6) {f,hoj, g )= Yl W/fo) - m)(g(v) - gt?)) + E y/C^K*) 



{x,y}CA L XGA L 
\x—y\=l 

for all f,g G £ 2 (A^). This shows that /iq.z, is non-negative. It also follows from (I2.5P and A.l 
that, uniformly in L, 

(2.7) \\h , L \\<4d\ max „ "^max 
Our second assumption is as follows. 

Assumption A. 2: We assume that /iq,l is positive definite. 

Since h 0t L > \ rmn xe ^ L k x , A. 2 holds for all L if there is a number k min with < fc mm < /c x . 
There are also other conditions which guarantee that A. 2 holds. Assumption A. 2 is sufficient 
to guarantee that Hl has purely discrete spectrum and a non-degenerate, normalized ground 
state Ql- This will be clear as a result of the diagonalization, which we continue with below, 
see e.g. Remark [2TT1 

The position and momentum operators q x and p x , on suitable domains, satisfy the self- 
adjoint form of the canonical commutation relations (CCR) 

(2.8) [Px,P y ] = [qx,q y ] = and [q x ,p y ] = i5 Xjy I. 

A canonical transformation is a linear change of position and momentum operators imple- 
mented via a symplectic matrix, see [18]. More precisely, define Q = (Qk), P = (Pk), 
fc = l,...,|Ai|,by 

(2.9) ( % ) = S - ( I 
where S is a 2|A/ / |-by-2|A/ / | symplectic matrix, i.e. 



(2.10) SJS 1 = J, J : 



-/ 

1 



One can check that the modified position and momentum operators still satisfy 
(2.11) [P 3 ,Pk} = [Q 3 ,Qk]=0 and [Qj, P k ] = iS j>k I, 
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and can be used to express the quadratic Hamiltonian in the form 
(2.12) H L = (Q T ,P T )S T U L s(^ Q p 

Every quadratic Hamiltonian of the form (12.31) with a real symmetric matrix Ml can be 
diagonalized in the sense that a symplectic matrix S exists such that S t MlS in (I2.12p becomes 
diagonal. This is the content of Williamson's theorem, e.g. Section 8.3 in [18]. In our case 
(12. 4p this can be done explicitly in two simple steps. 

For the first step let 

/i 1 / 2 
^' 2 



(2.13) Si 



One checks that Si is symplectic and that 



(2.14) s(m l s 



Tjrn C _ I h L 



/ 



where 

(2.15) h L = ^ 2 h ,LV 1/2 . 

By our assumptions, the real symmetric matrix hi is positive definite and satisfies the norm 
bound 

(2-16) \\h L \\ < IMHIMI < o ( 4cU 

max i ^"rnax/' 
Z77T. m ; n A 

In the second step /it, is diagonalized in terms of a real orthogonal matrix O, i.e., 
(2.17) T h L = 1 2 , 

with 7 2 = diag(7^) (and, for later reference, 7 = diag(7fc)). We have set 7^ to be the 
eigenvalues of hL counted with multiplicity, taking 7^ > for each k — 1, 2, . . . , |Ai|. 
Let 



O 
o 



(2.18) S 2 
which is symplectic. Thus 

(2.19) S = S 1 S 2 
is symplectic and 

(2.20) S T ( ho > L °)s-^ 2 ° 



^i 2 o 

/i" 1 / 2 o 



fi J \ I 

which we had set out to accomplish. 

In terms of Q and P from (12. 9p the harmonic Hamiltonian can be rewritten as 

(2.21) H L = (Q T ,P T ^ 12 



I I \ P 

(2-22) = £(P fc 2 + 7fc 2 QD- 
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This has the form of a system of uncoupled harmonic oscillators (or free Boson system) which 
can be explicitly diagonalized: Let 

(2.23) B k = x /?h.Q k + -± = p k , fe = l,...,|Ai], 

V 2 V 2 7fc 

a set of operators which together with their adjoints B k satisfy CCR in the creation-annihilation- 
operator form 

(2.24) [Bj, B k ] = [B*, B* k ] = 0, [Bj, B* k ] = SjjJ. 
One gets 

(2.25) H L = J2lk(2B* k B k + I). 

k 

It follows that 

• Hl has a unique normalized ground state Ql, with ground state energy 70 = Ylklk, 
which is characterized by B k Qi = for all k. 

• An orthonormal basis of eigenvectors for Hl is given by 

|A,| 

(2.26) j ( S ;)«*n L> 



• 1 V a i- 

3 =1 V J 

where atj, j G {1, . . . , |A^|}, are arbitrary non- negative integers. The corresponding 
eigenvalues are 

I A £ | 

(2.27) 7 o + 2^a j7i . 

i=i 

In particular, the ground state gap is 2minfc7fc. 

For many calculations, it is convenient to express the original position and momentum 
operators in terms of these operators B k which diagonalize Hl- Note that, with the explicit 
symplectic transformation in ( I2.19p . (12. 9p can be rewritten as 

q \ _ ( \(Q 

P J ~ V J\p 

Recalling (12. 23ft . we see that the S's are defined as 



(2.28) 



(2.29) 

and therefore 
(2.30) 



B \ J_f 7 1/2 iT 1 ' 2 \( Q 



q\ J_ / ^OT 112 \(\\\(B 

V ) ~ y/2 V /i" 1 / 2 ^ 172 ) V -i i ) V B * 

An explicit expression for the evolution of f!2.30j) under the dynamics generated by Hl will 
also be important. We briefly discuss this next. 

By Stone's Theorem, given an initial state ip e D(Hl), the unique solution of the Schrodinger 
equation 

(2.31) idtip{t) = H L ^(t) with ^(0) = ^, 
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is given by ip(t) = e~ ltHL tp, where the unitary operators e~ ltHh are defined via the functional 
calculus for self-adjoint operators. An alternative description of this dynamics is in terms of 
the so-called Heisenberg picture. Here, for example, one can define a one-parameter group of 
automorphisms r t on B(Hl) by setting 

(2.32) Tt (A) = e itHL Ae- ltHL for all A e B{U L ) ■ 

r t is called the Heisenberg dynamics generated by H^; it depends on L, but we suppress that 
in our notation. 

One can extend this notion of dynamics to certain unbounded operators on T-L^. In fact, 
appealing again to Stone's Theorem, one can check that the function bk(t) = Tt(Bk) is the 
unique solution of b' k {t) = —2ijkbk(t) with 6(0) = B^. (The previous equation holds, in the 
strong sense, on D(Hl) and the solution can be extended by linearity to the domain of B k .) 
In this case, we conclude 

(2.33) n(B k ) = e- 2i ^ kt B k and similarly T t {B* k ) = e 2i "< kt B* k . 
Applying fl2~33l) to (|2T30|) we find that 

Tt{q) \ J_f ^ l2 OT 112 \( e" 2 ^ e 2 ^ \ ( B 

nip) J ~ V2 V /i" 1/2 07 1/2 J V - %e ~ 2i ^ ie2tlt J\ B * 

We end this section with a few remarks. 



(2.34) 



Remark 2.1. While it is convenient to express the above reduction of Hl to a free Boson 
system in terms of the canonical transformation formalism, one may equivalently proceed by 
expressing the above transformations as unitary equivalences: Consider unitary operators U\ 
and U 2 on L 2 (IR Al ) defined by = cf{fi^ 2 q) and {U 2 f){q) = f(Oq). Here fi and O 

are the diagonal and orthogonal matrices defined above and c = Y[ x (^ m ^) 1 ^ normalizes U\ to 
be unitary. Let U := V^JJx- TJ\ rescales all masses to be 1/2 and U2 diagonalizes the quadratic 
potential, so that what we have shown above amounts to 

(2.35) UH L U* = -A + ^llll 

k 

meaning that Hl is unitarily equivalent to a system of uncoupled harmonic oscillators. This 
has eigenvalues given by with eigenvectors 112. 26}) taking the explicit form of Hermite 

functions, thus also showing completeness. From this point of view one gets the modified 
positions and momenta satisfying Ii2.22\) as 

(2.36) Q k = U*q k U, P k = U* Pk U. 

Remark 2.2. The above considerations can easily be modified to apply to more general qua- 
dratic Hamiltonians ( including multi- dimensional oscillators ) of the form 

hi \ / q 
h 2 \p 



T JT\ 



(2.37) H=(q 1 ,p 

with positive matrices hi and h 2 . Such cases are considered in [15]. The symplectic matrix in 
A2.19\) is generalized by 



(2.38) S 



h l 2 /2 O 
h~ 1/2 O 
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where O is the diagonalizing orthogonal matrix for h := h 2 h\h 2 ■ As above, the associated 
canonical transformation reduces H to a free Boson system. Furthermore, as in [15] or |16j . 
one could replace Z d by more general graphs and the Al by an exhausting sequence of finite 
subgraphs. 

We have chosen not to work in any of these more general settings, as our later applications 
to random oscillator systems will be stated in terms of the special case considered above. 
But generalizations of our results covering more general graphs or more general quadratic 
Hamiltonians are certainly possible. 

3. Dynamical Localization 

The first set of main results of our work is aimed at establishing a form of dynamical 
localization for the Heisenberg evolution f)2.32p of suitable observables under the Hamiltonian 
Hl. As discussed in the introduction, this will be done in the form of zero-velocity Lieb- 
Robinson bounds. 

While in this and the following sections we work with finite systems A^, the constants and 
decay rates in our bounds will be uniform in L. It is essentially for this reason that our proof 
of dynamical localization extends to infinite systems by taking the thermodynamic limit, as 
will be discussed in Section 

The mechanism leading to dynamical localization in oscillator systems will be disorder, 
introduced by assuming that some of the system parameters are random variables. Exponen- 
tial decay bounds will be proven for disorder averaged quantities. It is this disorder average 
which distinguishes our bounds from previously proven deterministic Lieb-Robinson bounds 
and allows to show that our bounds hold uniform in time, i.e. with group velocity zero. 

Our general strategy for studying the Heisenberg dynamics of Hl is to reduce it to the 
dynamics of the operator Hl defined by f)2.15p . which will serve as an effective one-particle 
Hamiltonian for the many-body Hamiltonian Hl- We will state our main results, in this 
section as well as in Sections S] and as saying that localization properties of Hl imply 
localization properties of Hl- In Section [H] we will provide regimes in which the required 
localization properties of Hl can be verified. 

In this section we will study the Heisenberg dynamics for two types of observables. Most 
natural may be to consider the dynamics of the position and momentum operators q x and p x . 
We will postpone this to Section I3.2[ while we will first study the dynamics of Weyl operators 
associated with positions and momenta. The main reason for this is that the Weyl operators 
are bounded (in fact, unitary) which will lead to a slightly more streamlined presentation of 
our results. 

Throughout this section as well as in Sections H] and we will assume A.l and A. 2 without 
further reference. 

3.1. Dynamical Localization for Weyl Operators. In this section, we introduce the set 
of Weyl operators W(f) generated by functions / : — > C. In terms of these observables, 
we prove a strong form of dynamical localization for certain random oscillator models, see 
Theorem 13.31 and the comments that follow. 

To define Weyl operators, it is convenient to introduce annihilation and creation operators. 
These are defined, for any x G Ai, by setting 

(3.1) a x = —= (q x + ip x ) and a* = —= (q x - ip x ) . 
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Recalling (I2.8p . it is easy to see that these operators satisfy the CCR, i.e., 

(3.2) [a x ,a y ] = [a*, a*] = and [a XJ a* y ] = 5 X:V I . 
More generally, for any / : — > C, one can set 

(3.3) a(/) = ^ Jga, and a*(f) = 

It follows from (13. 2 j) that 

(3.4) [a(/),a(y)] = [fl*(/),a*(y)] = and [«(/), «*(<?)] = {f,g) . 
Each /, as above, corresponds to a Weyl operator W(f) defined by 



(3.5) 



W{f) = exp 



^ («(/) + «*(/)) 



exp 



(Re[/(x)]g x + Im[/(x)]p a 



where the second identity is (13. 9 p below. Clearly, is a unitary operator on T-Ll with 

= W(— /) and W(0) = I. As is well-known, see e.g. [10], Weyl operators satisfy the 
so-called Weyl relations, i.e. for any such / and g 

(3.6) W(f + g) = e-2 Im ^'^W(f) W(g) . 

These Weyl operators depend implicitly on the volume Al, but we suppress this in our nota- 
tion. 

A convenient deterministic fact about these Weyl operators is that they set, 
invariant under the harmonic dynamics. For completeness, we state this as a lemma. 

Lemma 3.1. For each finite set Al, let Tt denote the corresponding harmonic dynamics. For 
any f : A L C and t6K, 

(3.7) Tt (W(f)) = W(f t ) where f t = V^e^V f 
with a real linear mapping V defined by 

(3.8) Vf = 7 - 1/2 OV /2 Re[/] + * 7 1/2 0V 1/2 Ihi[/] , 
and 7, O, and fi are as in Section^ 

Proof. To prove this fact, we first express the Weyl operators in terms of the operators Bj, that 
diagonalize the Hamiltonian, see (I2.25p . Using vector notation, write a = (a x ) and a* = (a*) 
and observe that 



(3.9) 



a(f) + a*(f) = (J J 1 ) 



a 

a* 



1 % 
1 -% 



Appealing to (12. 30 p . it follows that 
(3.10) 

a(f) + a*(f) = (Re[/] T ,Im[/n 



((Vff, (Vf) 



V2(Re[/] T ,Im[/] 



fjl/2 0l -l/2 o 

o /i-^oy/s 

B* )=B(Vf) + B*(Vf) 



t\ ( q 
v 



l l 

— i i 



B 

B* 
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with a real-linear mapping V given by (I3.8p . Here, in analogy to (13. 3p . we have also introduced 
the notation 

(3.11) B^) = YJtF)B k and B\<j>) = ^ 0(fc) B* k . 

k k 

We have shown that 

i , , . *, ,u' 

: exp 



(3.12) W(f) = exp 



(a(f) + a*(/)) 



The claim in (13. 7p now follows from (I2.33p . □ 

Before we state our first result, the following observations will be useful. It is easy to check 
that the inverse of the mapping in (13.81) is given by 

(3.13) V~ x <\) = /i" 1/2 C 7 1/2 Re[0] + i/i 1 / 2 O7- 1 / 2 Im[0] , 
in which case, one readily verifies that 

(3.14) Re[/ t ] = /i" 1/2 cos(2t^ /2 )/i 1/2 Re[/] - ^ 1,2 h\ /2 sm(2th L /2 )fi- 1/2 lm[f] 
whereas 

(3. 15) Im[/ t ] = /i 1/2 cos(2t^ /2 ) y u- 1 / 2 Im[/] + tL 1/2 h~ L 1/2 sm(2th L /2 ) /iV 2 Re[f] . 

Here, and below, functions of the operator Hl are defined via the elementary functional calculus 
of self-adjoint operators in finite-dimensional Hilbert spaces, e.g. cos(2t/i^ 2 ) = Ocos(2t7)0 T , 
and so on. Note that, due to (I2.17p . we have O^O t = h 1 / 2 , which is the reason for the 

1/2 

frequent appearance of h L in relations between the one-particle operator hi and the many- 

1 12 

body operator Hl, here as well as below. In fact, it is more accurate to think of h£ rather 
than of Hl as the effective one-particle Hamiltonian associated with H^. 

The identity (13.71) in conjunction with ( 13. 14ft and ( 13.151) allows us to derive dynamical 
localization properties of Hl from localization properties of Hl- All our considerations so far 
have been deterministic, meaning that they hold for any choice of the parameters k x , m x and 
\ XtV satisfying our basic assumptions A.l and A. 2. For our results on localization we will 
think of these parameters as random variables on a probability space (Q, P), which satisfy A.l 
and A. 2 almost surely. By E(X) we denote the expectation (average) of a random variable 
lonfi with respect to P. 

For the sake of definiteness, one may think of the special case where one or several of the 
three parameter sets are given by i.i.d. random variables. In fact, our applications in Section [6] 
will be for the case where the k x are i.i.d., while the other parameters will be kept deterministic 
(for example, constant). But our general results below would, at least in principle, allow for 
applications with other types of disorder. 

The required one-particle localization property will generally be expressed in the following 
form: 

Definition 3.2. We say that Hl has exponentially decaying (a,r)-eigenfunction correlators if 
there exist C < oo and \i > such that 

(3.16) E ( sup \ (5 x ,h a L u(h)5 y \ r ) < Ce~^ x - y \ 

\H<i / 
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for all positive integers L and all x,y G A^. Here the supremum is taken over all functions 
u : R — > C which satisfy the pointwise bound \u\ < 1. 

As discussed in the introduction and, more thoroughly in Appendix |A] eigenfunction cor- 
relators for a = are a well-known object in the literature on Anderson localization. A new 
aspect of our work is that we will need them for other, in particular negative, values of a. 

Our first result for disordered oscillator systems establishes dynamical localization for the 
Heisenberg evolution of Weyl operators. 

Theorem 3.3. Assume that hi has exponentially decaying (— 1/2, r)- eigenfunction correlators 
for some r G (0, 1]. Then there exist C and /i' for which 

(3.17) e( sup \\[T t (W(f)),W(g)}\\) <C , ^\f(x)\ r \g(y)\ r e-'<M 
for any L > 1 and f, g : Ai — > C. 

Proof. Let / and g be as above. Using Lemma [3.11 and the Weyl relations ( 13. 6ft . it is easy to 
see that 

(3.18) [t*(W(/)), W(g)} = ( e -*M</^>] - l) Wfo)W(/ t ) . 
In this case, the norm bound 

(3.19) l|[n(W(/)),W(p)]|| <min{2,|Im[(/ f ,4>]|} 
readily follows since the Weyl operators are unitary. Moreover, since 

(3.20) M(f t ,g)] = (Re[f t },Im[9}) ~ (lm[f t ],Re[g}) , 
a short calculation, using (I3.14p and A3. 15j) . shows that 

(3.21) 

lm[(ft,9)\ 

= (^ 2 Re[f],hl 1/2 sm(2th]! 2 )^ 2 Re[g}) + (^ 2 lm[f], cos(2< /2 y/ 2 Re[<7]> 

+ ( / i 1 / 2 Re[/],cos(2< /2 ) / i- 1 / 2 Im[ 5 ]) - (^ 1/2 lm[fl hf sin(2< /2 )/i" 1/2 ReM) . 

By A.l the entries of the diagonal matrices /x 1//2 and /x -1 / 2 satisfy almost sure L-independent 
constant bounds. As a consequence, almost surely, 

(3.22) 

S u V \\[r t (W(f)),W(g)}\\ < sup |Im[(/„y>] | r 

< c^i/^n^r E su P i(4,^(wr- 

a:,?/ ae{-l/2,0,l/2} 

The supremum over t G 1 has been absorbed into the supremum over |u| < 1. Among the 
three values of a on the right of ( I3.22p . the term with a = —1/2 is the most singular and 
crucial one. In fact, if |u| < 1 and u(s) := s 1 ^ 2 u(s), then \u\ < (max{<r(/iL)}) 1//2 on the 
spectrum a (hi) of hi. From this it readily follows that 

— 1/2 

(3.23) sup \ (5 x ,u(h L )5 y )\ < C 2 sup \ {S x ,h L u(h)5 y )\, 

\u\<l \u\<l 



OSCILLATOR SYSTEMS 



13 



where C2 can be chosen as the square root of the constant on the right hand side of (12.16p . 
In a similar way, the a = 1/2 term in (13.221) can also be absorbed in the a = — 1/2 term. The 
claim in ( 13 . 1 Tf) now follows from ( 13 . 1 6f) with a = — 1/2. □ 

Two remarks are in order to clarify the relation of our assumptions to the different a-terms 
in ( m> : 

Remark 3.4. As is clear from \3.2lS) , the most singular term a = —1/2 in Ii3. 22\) is due to 
the real parts of f and g, which, by / Iff. 5\) and 113. 9) . correspond to contributions by positions 
to the Weyl operators. If f and/ or g are purely imaginary (corresponding to Weyl operators 
only involving momenta), then a bound of the form \3.17) can be proven under (for the sake 
of applications) weaker assumptions than exponential decay of (—1/2, r) correlators. For ex- 
ample, exponentially decaying (1 / 2, r)-eigenf unction correlators for some r e (0,1] suffice to 
conclude that 

(3.24) E f sup||[Ti(e^),e^]|| j < CV' 1 ^ 1 . 

This non-symmetry of our results in the positions and momenta is, of course, due to our 
initial choice of the model \2. where interactions are introduced via the position operators. 

Remark 3.5. That in the above proof the a = and a = 1/2 terms could be absorbed into 
the a = —1/2 term was due to the fact that our assumptions give a norm bound on which 
is uniform in the disorder as well as in the size of the system, thus avoiding that the a = 1/2 
term introduces a singularity at infinite energy. Similarly, we could avoid the zero-energy 
singularity arising from the a = —1/2 term by, in addition to A.l, requiring a lower bound 
k x > ^min > ; uniformly in x and the disorder. This would show hi > k m i n /2 and, via Ii2.15\) 
and \2.27\) , result in a uniform ground state gap for the oscillator system Hl (with an explicit 
lower bound given by (k m - m / , m mffl ) 1 '' 2 ) . In this case all (a,r) correlators would be equivalent 
in the sense that their exponential decay would not depend on the choice of a. 

The following corollary is an immediate consequence of (13. 17j) in cases where / and g have 
disjoint support. It phrases our result in a form more similar to traditional Lieb-Robinson 
bounds and explains why we have referred to it in the introduction as a zero-velocity Lieb- 
Robinson bound. 

Corollary 3.6. Fix two finite sets X, Y C Z d with X n Y = 0. Let L > 1 be sufficiently 
large so that X U Y C A^. If hi has exponentially decaying (— 1/2, r)- correlators for some 
r G (0, 1], then there are C < 00 and fi' > such that 

(3.25) E(snp\\[r t (W(f)),W(g)]\\] < C'\\fUgle~^ x ^ 



holds for all f and g supported in X and Y respectively. Here d(X, Y) = mm x£ x,y&Y \x — y\. 

In [33] oscillator systems of the form ( 12. ip with constant parameters m x = m, 
^x,x+e- — ^j, j — l,---,d, were considered and, in particular, a deterministic LR-bound of 
the form 

(3.26) ||[T t (W(/)),W/(s)]||<C||/|U^ 

was proven for this case, see Corollary 3.2 in [33]- The finite positive number v is interpreted as 
a bound on the velocity of propagation of signals in the oscillator system. In the cases discussed 
in Section [6] where we can verify exponential decay of (— 1/2, r)-eigenfunction correlators for 
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Hl and after disorder averaging, our result shows that this velocity is zero for the random 
oscillator systems considered here. 

More general deterministic oscillator systems (allowing for variable coefficients and more 
general lattices) were subsequently considered in [IE] , again establishing a "finite speed of 
light" in the system. We stress here that there were no previous results establishing dynamical 
localization of oscillator systems in the sense of a zero-velocity Lieb-Robinson bound such as 
(153511 . 

In addition to the dynamics of Weyl operators, the papers [33] and [16] also have results for 
the dynamics of positions and momenta. The effects of disorder on the latter are the topic of 
the next subsection. 

3.2. Dynamical Localization for Positions and Momenta. Theorem 13.31 was concerned 
with the dynamics of Weyl operators. We will now consider an analogue for the dynamics of 
position and momentum operators. Due to the unboundedness of these operators, the results 
are a bit different. We comment on these differences in this section. 

We begin with a well-known result which we state as a lemma and prove for completeness. 
For any x,y G and t G R set 

[n(qx),q y ] [n(qx),P y ] 
Vt{Px)Ay\ [n{Px),Py] 

where t% is the harmonic dynamics, see (12.321) above. Our normalization is such that 

/ 



(3.27) A x , y {t) 



(3.28) A Xty (0) = 5 : 



.i.y 



-I 



Lemma 3.7. For any x,y G and each t G R, one has that 
(3.29) 

, ,., ( -(/i 1/2 5 a; ,/iZ 1/2 sin(2< /2 ) / i 1 /2^)/ ( / uV2^ )COs(2 t/ i V 2 ) / ,-i/2^)/ 
x ' y[ 1 \ -{^ l lH x ,cos{2th^ 2 )^H y )I -( /U - 1 /2^ ; / l V 2 siri (2t4/ 2 )^-V2^) / 

where /i and are as in Section 

Proof. The commutators appearing in A Xiy (t) are entries in the larger matrix 

(3.30) -((^)(r,fef,n(pn) T 

Note that the CCR for the -B^'s can be expressed as 

(3.31) ( I ) (B T , (B") T ) " (( |. ) (B T ,(B-f))= ( °j o ) • 

Using (13.311) . a short calculation shows that the matrix in (I3.30P simplifies to 
(3.32) 

^/2 Ql -x/2 q \ / sin (2 7 t) -cos(2 7 t) \ / 7 - 1 / 2 T /i 1 / 2 

^O^' 2 J V cos(2 7 t) sin(2 7 t) J \ 1 x I 2 O t ^ 

fi^O-y- 1 sin(2 7 t)OV /2 ~V 1/2 cos(2 7 t)G T /i" 1 / 2 
yT^O cos(2 7 t)(9 T /i 1 /2 M -V2 07 sin(2 7 t)0 VV2 

The claim in (13.291) now follows by taking the appropriate matrix entries. □ 
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In applying Lemma 13.71 to disordered oscillator systems, unlike with the Weyl operators, 
the upper left entries of (I3.29p . corresponding to the position-position case, do not satisfy a 
norm bound which is uniform in the system size and the disorder (and thus we can not use 

________ 1/2 

an argument as in the first line of ( I3.22p ). This is due to the appearance of the factor h L , 
which can have arbitrarily large norm as hL may have eigenvalues arbitrarily close to zero. It 
is for this reason that the following result is stated in two parts, with an additional fractional 
power included in the commutator bound for the position-position case. Its proof follows 
immediately from Lemma 13.71 For the sake of simplicity we do not state a third separate 
result for the momentum-momentum case, where the assumption could be weakened further 
by assuming exponential decay of (1/2, l)-eigenfunction correlators. 

Theorem 3.8. (a) If hL has exponentially decaying (— 1/2, r)-eigenj 'unction correlators for 
some r 6 (0, 1], then there exist C < oo and fi' > for which 

(3.33) E (sup \\[r t (q x ),q y }\\ r ) < C'e^'^ 

for all L and i,?/GAj,. 

(b) If hL has exponentially decaying (0, 1) -eigenf unction correlators, then there exist C < oo 
and // > for which 

(3.34) E fsupmax{||h(p x .),g,]||, \\[r t (q x ),p y }\\, \\[r t (p x ) , p y ]\\}) < C'e^'^ 

for all L and x, y G Al- 

We conclude this section by noting that it is tempting to try to derive Theorem 13.81 as a 
consequence of Theorem l3.3l at least if the assumption on eigenfunction correlators holds with 
r = 1. Considering only position operators, (13.171) would then imply that 

for every e > 0. Formally, this yields (I3.33P as e — > 0. However, as the position operators are 
unbounded, it is not obvious how to rigorously justify this limit. 

4. Ground State Correlations 

Localization near the bottom of the spectrum is sometimes described as the existence of a 
"mobility gap" , which does not mean that there is an actual gap in the spectrum above the 
ground state. In this section we show that in the models we consider a mobility gap implies 
exponential decay of spatial correlations in the ground state. 

As is discussed in Section [2] after (I2.25p . the oscillator hamiltonians Hl each have a unique, 
normalized ground state VLl which is characterized by B^l — for all k — 1, 2, . . . , |A^|. To 
fix notation, for any linear operator A in Hl such that Ql £ D(A), the domain of A, let us 
denote by 

(4.i) (A) = (n L ,An L ) 

the ground state expectation of A. In this section, we will investigate static ground state cor- 
relations (AB) — (A) (B) as well as dynamic ground state correlations (r t (A)B) — {r t (A))(B) = 
(r t {A)B) — (A)(B). As in the previous section, we have results for the case where A and B 
are Weyl operators as well as for the case of positions and momenta. 



(3.35) 









I sup 













-,ieq y 
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4.1. Correlations of Weyl Operators. We start with ground state correlations of Weyl 
operators. As is well-known (or follows by taking the /? — > oo limit in (I5.2p below), for any 
/ : Al C 

(4.2) (W(f)) = e-^ll^ll 2 

where V is the real-linear mapping given by ( 13. 8p . Let us denote the dynamic ground state 
correlations of Weyl operators W(f) and W(g) by 

(4.3) C(f, g; t) = (r t (W(f))W(g)} - (W(f))(W(g)) . 

In this context, our result is the following. 

Theorem 4.1. If has exponentially decaying (— 1/2, r)-eigenf unction correlators for some 
r G (0, 1], then there exist C < oo and ji' > such that 

(4.4) mLn V \C{f, g] t)\) <C'J2\f(x)\ r \9(y)\ r e-^ x - yl 
/or any f,g E f(A L )- 

Proof. Using Lemma [3.11 and the Weyl relations (13. 6p . one sees that 

(4.5) (r t (W(f))W(g)) = (W(f t )W(g)) 

= e-^ m[{ft ' 9)] (W(f t + g)} 
= e -|im[(/ t , 9 >] e -i||v(/ t +3)|| 2 

where we also used ( 14. 2p . Moreover, one has that 

(4-6) \\V(f t + g)\\ 2 = \\Vf t \\ 2 + \\Vg\\ 2 + 2Re[<V/ t) Vg)} . 

We have shown 

(4.7) C{f, g; t) = ^ e -^^)] e ^M(vfuV 9 )] _ A e -Kl|v/|| a +l|v fl ||") 

Now for any real a and b 

(4.8) e m+b - 1 = e\e ia - 1) + e b - 1 and so \e m+b - 1| < (\a\ + |6|)e |b| . 
In this case, 

(4.9) \C(f,g;t)\ < \ (|Im[</*, g}}\ + \Re[(Vf t , Vg)}\) ell^(^^)]l e -Kll^ll 2 +ll^ll 2 ) . 
Since it is also true that 

(4.10) \Re[(Vf t ,Vg)]\ < \\Vf t \\\\Vg\\ = \\Vf\\\\Vg\\ < \ {\V ff + \\Vg\\ 2 ) , 
the bound 

(4.H) W,g;t)\ < \ (|Im[(/ f ,o}]| + \Re[(Vf t ,Vg))\) 

follows. 

The first term above is identical to the term appearing in ( 13.191) . and we estimate it as 
before by using (I3.2ip . 
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The second term involves 

(4.12) 

Re[(Vf t ,Vg)} 
= (Re[e 2i *Vf], Re[Vg]) + (lm[e 2 ^V f},lm[Vg]) 
= (cos(2 T t)Re[V7], MVg}) ~ (sm{2jt)Im[Vf), Re[Vg}) 

+ (sm(2^t)Re[Vf],Im[Vg]) + (cos(2~/t)lm[V f], lm[Vg]) 
= (^ 2 Re[f], hi 1 ' 2 cos(2< 2 ) / i 1 / 2 ReM) - (fi~ 1/2 lm[f], sin(2< /2 y^Re^]) 
+ ( / i 1 / 2 Re[/],sin(2t4 /2 ) / i- 1 / 2 Im[^]) + (/i^M/] , /if co^th 1 ' 2 )^ 1 ' 2 !^]) 

which is very similar to ( I3.2ip . the main distinction being that the role of sines and cosines is 
reversed. As \C(f,g;t)\ is uniformly bounded by 2, we can argue as in H3.22j) to introduce a 
fractional power r G (0, 1]. As in the proof of Theorem 13.31 the bound (I4.4p now follows from 
the assumption on decay of eigenfunction correlators. 

□ 

For the averaged static correlations of Weyl operators, e.g. t = in (14. 3p . one extracts the 
following explicit bound from the proof of Theorem 14.11 

Corollary 4.2. The static correlations of Weyl operators satisfy the bound 

(4.13) E(\C(f,g;0)\) 

< l\lm[{f,9)]\+C^2\f(x)\\g(y)\ (e(K*,,/#\)|) +E(\(5*,hi 1 \)\j) 

for all L and f, g : A^ — > C. Here C = | max{m max , 1/^mm} for the constants from A.l. 

This will allow to get exponential decay bounds on static correlations for a wider range 
of applications than on dynamic correlations. In particular, exponential decay of the terms 
K(\(S X , h 1 / 2 S y )\) and ~E(\(5 X , h L 1 ^ 2 5 y )\) will only require that the effective one-particle Hamil- 
tonian is localized near energy E = 0, while exponential decay of eigenfunction correlators 
( I3.16P requires localization of hi at all energies, see Section [6] and Appendix lAl 

4.2. Correlations of Positions and Momenta. We next prove localization-type results for 
the ground state correlations associated with the position and momentum operators. 

Let Oj(x) be the matrix elements of the orthogonal matrix O from (12.1 7p with respect to 
the canonical bases in £ 2 (Al) and £ 2 ({1, . . . , |A^|}). Then the components of (12.301) are given 
by 

(4-14) q x = \m-V 2 £ O k (x) lk ~ 1/2 (B k + BJ) 

k 

and 

(4.15) p x = ivr^J 2 O k (xhl /2 {-B k + Bl) . 

k 

for each x G A^. Using these relations and the fact that B^l = for all k, we find that 



(4.16) 



(Px) = (q x ) = for all x e A L 
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For the Heisenberg evolution of the positions and momenta we have from (12.341) . 

(4-17) r t (q x ) = \mZ 112 £ O k (x)^ 1/2 (e^B k + e 2 ^B* k ) 

k 

and 

(4.18) r t (p x ) = im l J 2 O k {xW,J 2 {-e^B k + e»»«SJ) . 

k 

For given x and y in A^,, using (I4.16p . we can express the dynamic ground state correlations 
of positions and momenta at x and y, respectively, in terms of the effective one-particle 
Hamiltonian hi as follows: 

Lemma 4.3. For any x, y G and t G M., 

(4.19) (r t {q x )q y ) = ^Km^ 1 /^, h~ L 1/2 e~ mh ^ 8 y ) , 

(n(q x )p y ) = '-(m y /m x )y 2 {S x ,e- 2 < /2 S y ), 

(n('Px)q y ) = ~(m x /m y ) 1/2 (5 x ,e- 2ith L 2 5 y ) , 

(r t (p«)p y > = M^ii/ife^ 2 ^. 

Proof. We calculate the first identity. The rest are done similarly. Combining ( 14.14p with 
(I4.17P and using that the B k annihilate the ground state and that B^VLl is normalized by 
( I2~26j) . one gets 



(4.20) (r t (q x )q y ) = -(m^r^^k^OM^ 2 ^ 

k,k' 



1/2 -1/2 

X 



x ((e- 2i ^B k + e^Bl) (B k , + B* k ,)) 
\{m x m v )-^ £ O^O^yh-^-^e-^iB.B;,} 



k,k 

1 



1 / 2 ^O fc (x)O fc ( 2 /) 7 ,: 1 e- 2 ^' 



it 

= ^(m.rn,)- 1 / 2 ^, Oy-ie-'VCFSy), 
which gives (j4.19|) . □ 
An immediate consequence is 

Theorem 4.4. (a) If hz has exponentially decaying (— 1/2, r)-eigenj 'unction correlators for 
some r G (0, 1], then there exist C < oo and // > stzc/i i/iai 

(4.21) E (sup \{n(q x )q y )\ r ) < C'e'^^ 

\tett J 

for all L and x, y G A^. 
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(b) If h L has exponentially decaying (0, l)-eigenfunction correlators, then there exist C < oo 
and fjf > such that 

(4.22) E (supma^KT^)^}!, \(rt(q x )p y )\,\(r t (p x )p y )\}] < C'e^'^ 



for all L and x,y G A^. 

As for the Weyl operators in Section 14.14 we a ls° get simple explicit bounds for the static 
correlations of positions and momenta: 

Corollary 4.5. We have (q x p y } = —(PxQy) — \$x,y as well as 

(4.23) E(|(g^)|) < ^E(|(^,^ 1/2 ^)|) 
and 

(4.24) m<PxP v )\) < m^Ed^,^ 2 ^)!). 

5. Correlations for Thermal States 

On physical grounds (increasing entropy) one would expect the correlation length to de- 
crease with increasing temperature. We therefore expect that the conditions under which we 
can prove exponential decay in the ground state, would also be sufficient to prove exponential 
decay at positive temperatures. This is not entirely trivial however. In Section [3] we proved 
exponential decayin the ground state based on localization near the bottom of the spectrum. 
At higher temperatures higher portions in the spectrum dominate the state and the corre- 
sponding states may be delocalized. It is therefore not entirely obvious that exponential decay 
holds for all temperatures. Nevertheless, in this section we prove exponential decay of spatial 
correlations in the equilibrium states (also called thermal states) at any temperature. 

Thermal states are defined as follows. Fix e (0, oo). For any A such that Ae~ f5HL is trace 
class, set 

Tr \Ae~ pHL ] 

to be the expected value of A in a thermal state corresponding to a positive temperature 
inversely proportional to /3 < oo. Note that, as for zero temperature, (jt(A))p = (A) p. 

5.1. Correlations of Weyl Operators. We start with the calculation of thermal state cor- 
relations of Weyl operators, which provides a positive temperature analogue of the results in 
Section O 

It is well-known, see Proposition 5.2.28 of |10] or Chapter XII. 12 of [31J . that 

(5.2) {W{f))f, = e-t\\ Al '" v, ( 
where V is as in (I3.8P and 

(5.3) A p = (I + e- 2 ^)(I - e" 2 ^)" 1 = coth(/3 7 ). 

Here 7 is the diagonal operator from (I2.17p . 

The thermal correlation of the Weyl operators W(f) and W(g) is 

(5-4) C fi (J, g; t) = (r t {W{f))W{g)) p - (W(f))p(W(g))p . 
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Moreover, the Weyl relations (13 .6p . (13 .7p and (15. 2p imply that 

(5.5) (r t (W{f))W{g)) p = e -i Im ^^] e -ill<V(/ t+9 )|| 2 _ 
Like before, the relation 

(5.6) \\AfV{f t + g)\\ 2 = \\AfVf t \\ 2 + \\AfVgf + 2Re[{A f) Vf t , Vg)\ 
holds and therefore we obtain 

(5-7) n(f,g;t)\ < \ (|Im[(/ t , g))\ + \Re[(A p Vf t , Vg)}\) 

as in Section S7JJ By a calculation generalizing (I4.12j) we get 
(5.8) Re[(^V/ tj ^)] 

= (^ 2 MflMhL)» 1/2 Mg\) - {^MflfoM^Reig}) 

+ (^ 2 Mflfo,t{hL)»- ll2 Mg\) + (M- 1/2 Im[/],0 3 ,*(^)/i- 1/2 ImM). 

Here we have set 
(5.9) 

M ( S ) = s ~ 1/2 coth(/3s 1/2 ) cos(2ts 1/2 ), <j) 2tt (s) = coth(/3s 1/2 ) sin(2ts 1/2 ), 3jt (s) = sfa^s). 

The operator functions 0i,t(/iL), <t>2,t{h>L) and appearing here are more singular than 

the operator functions we have encountered earlier in (13.211) and (]4.12p . In fact, due to the 
additional factor coth(/3s 1 ^ 2 ) which behaves like l/(/3s 1 ^ 2 ) near s = 0, we have 4>i,t(s) ~ l/(fis) 
for the most singular function in (15.91) . 

The stronger singularities are the mathematically most interesting new feature appearing 
in our treatment of thermal state correlations, as opposed to the earlier results on Lieb- 
Robinson bounds and ground state correlations. In the applications in Section |6] below we 
will have to control more singular types of eigenfunction correlators of the effective one-particle 
Hamiltonian hi than what is needed in applications of our earlier results. 

Under the assumption that the required more singular eigenfunction correlators are expo- 
nentially decaying, we thus get the following analogue of Theorem 14.11 for the thermal state 
correlations of Weyl operators. 

Theorem 5.1. If Kl has exponentially decaying (— l,r)- eigenfunction correlators for some 
r G (0, 1], then there exist C < oo and fi' > such that 

(5.10) E(snp\C p (f,g;t)\] < C'J^ \f(x)\ r \g(y)\ r e^^ 

for any f,g E i 2 (A L )- 

Remark 5.2. In principle, one can derive Theorem 4-1 from Theorem I5.il by taking the 
(3 — > oo limit. This would, however, require us to work under the stronger assumption of 
(—1 , r) - correlator decay instead of (—1/2, r)- correlator decay as in Theorem J^.l, which seems 
unnatural and weakens applicability. For this reason, we have separated our results into two 
sections. 



As in Section 14.11 we state an explicit bound for the static thermal correlations of Weyl 
operators, which will prove their exponential decay in a wider range of applications. 



OSCILLATOR SYSTEMS 21 



Corollary 5.3. The static correlations of Weyl operators satisfy the bound 
(5.11) 

EflcjK/^o)! 1 /*) 



< 

- 2V2 



1/2% 



+ E(\(5 x ,h~ L 1/2 coth^f)^)! 1 / 2 ) 

/or a// L and /, g : Ax — > C. iJere C = |max{m max , l/m min } /or i/ie constants from A.l. 

5.2. Correlations of p's and g's. Finally, we consider thermal correlations of positions and 
momenta. For them we have the following deterministic facts. 

Lemma 5.4. Fix f3 > 0. One has that 

(5.12) (q x )i3 = ij) x ) p = for all x G A L . 
Moreover, for any x, y G A^ and t el, 

(5.13) (T t (qx)q y )/3 = ^{m x m y )- 1/2 (5 x ,^(p lt t(h L ) -ihl 1/2 sin(2th^ 2 )^d y ) 
(n(qx)Py)/3 = ^(m y /m x ) 1/2 (5 x , jcos(2t/^ /2 ) - i(h,t( h L)j S v ) 
(n(Px)q y )p = -^(m x /m y ) 1/2 (5 x , |cos(2t/iJ / /2 ) - i(j) 2 ,t{h L )} S y ) 

(n(Px)p y )p = (m x m y ) 1/2 (5 x , \j>u{h L ) - ih 1 / 2 sin(2t/ij / /2 )| S v ) , 
where we use the functions defined in $5.9\) . 



Proof. This is proven by calculations similar to those used to prove Lemma l4~3l this time using 
the following facts on the thermal state expectations of the bosonic creation and annihilation 
operators: 

(5.14) (B k ) p = {Bl)p = {B k B k ,) p = (BtB* k ,) = , 

(5.15) (B k B* k ,) p = (1 - e- 2 ^)- l 8 k>w , (B* k B k ,) p = (e 2 ^ - 1)'%^ , 
see Proposition 5.2.28 of [TO] . 

□ 

_______ _______ 1/2 

The crucial difference between ( 14.1 9p and (I5.13P is that in the latter the factor coth(/3/i L ) 
in 0i,j(/il), (p2,t(^L) and <t>z,t{hLj leads to a stronger zero-energy singularity, in each case 

— 1/2 

corresponding to an extra factor h L . Thus in each of the two parts of the following thermal 
state analogue of Theorem 14.41 a correspondingly more singular eigenfunction correlator needs 
to be controlled in the assumption. Otherwise the claims follow immediately from Lemma l5~4l 

Theorem 5.5. (a) Ifh^ has exponentially decaying (—1, r)- eigenfunction correlators for some 
r G (0, 1], then there exist C < oo and fi' > such that 



(5.16) 

for all L and x, y G A^. 



E(sav\{r t (q x )q v )pA <C'e^'^ 
\tm J 
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(b) If Hl has exponentially decaying (—1/2, r)-eigenf unction correlators for some r G (0, 1], 
then there exist C < oo and yi! > such that 

(5.17) E ( B n V max{\(T t {p x )q y ) p \ r , \(n(q x )p y )p\ r , \(r t (p x )p y }^}) < ^'^ 

for all L and x, y G Al- 

As before, for the static correlations we state an explicit bound tailored to our applications 
in Section [HJ 

Corollary 5.6. We have {q x p y )i3 = —(Px<ly)i3 = \$x,y as as 

(5.18) E(\(q x q y )^ 2 ) < -^E(\(6 X , Z^ 2 coth^f )5 y )\ 1 ' 2 ) 

and 

(5.19) E(\{p xPy )p\) < m max E(|(^,4 /2 coth(^ /2 )5,)|). 

6. Applications 

Our results on disordered oscillator systems Hl in Sections El H]and|5]were stated in terms of 
assumptions on the exponential decay of eigenfunction correlators (13 . 161) of the associated one- 
particle Hamiltonian hi. Thus applications of our results will consist in providing conditions 
on the choice of the random parameters m x , k x and \ XtV under which decay of eigenfunction 
correlators can be verified. 

The main application will be what we will refer to as the Anderson case, which will describe 
the following set of conditions: 

Assumption A. 3: Let m x — m G (0, oo) and X x>y = A G (0, oo) be constant, and k x , 
x G Z d , be i.i.d. random variables with absolutely continuous distribution given by a bounded 
density p supported in [0, /c max ], < A; max < oo. 

With this choice, the operator Hl in (12.151) . up to constant rescaling, becomes the Anderson 
model discussed in Appendix |A] This makes a wealth of known localization properties avail- 
able. In particular, exponential decay of (0, l)-eigenfunction correlators (see Definition 13. 2p is 
well known in various regimes. In fact, it is in works on the Anderson model where the term 
eigenfunction correlator was first used in this context, e.g. [I], [1] and also the recent review 

52]. 

However, singular eigenfunction correlators (a < in (I3.16P ) do not seem to have appeared 
in previous works. Establishing localization bounds for them, after reviewing known results, 
is the main content of Appendix |A] below. What makes the extension of localization results 
for eigenfunction correlators to singular eigenfunction correlators non-trivial is the fact that 
in A. 3 we allow the k x to take values arbitrarily close to zero. As discussed in Remark 13.51 

— 1/2 

above, this means that the operator h L is not uniformly bounded in the system size and 
the disorder and thus eigenfunction correlators for different values of a are not equivalent. 

Expressed in terms of the many-body Hamiltonian Hl this means that in cases where ex- 
ponential decay of eigenfunction correlators can be shown under Assumption A. 3, we get 
many-body localization properties (e.g. zero-velocity LR-bounds as well as exponential clus- 
tering for ground states and thermal states) in situations where Hl does not satisfy a stable 
ground state gap condition. 
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We group our applications to the Anderson case into three theorems, corresponding to the 
three different regimes described in Proposition IA.1I where localization is known. We start 
with the general case A. 3, without restriction of the dimension, where localization is known 
to hold near the bottom of the spectrum. This does not allow to prove exponential decay of 
eigenfunction correlators, as required by all of our results in Sections [31 |4] and El that refer to 
dynamical quantities. However, localization of Hl near the bottom of the spectrum suffices to 
show exponential decay of the static (t = 0) ground state and thermal state correlations for 
Weyl operators as well as positions and momenta: 

Theorem 6.1. Assume A. 3 and fix (3 e (0, oo). Then there exist constants C < oo and fi > 
such that 

(6.1) E(\C(f,g;0)\)<Cj2\f^)\\9(y)\e-^ x - yl 
and 

(6.2) E(\C p (J,g;0)\ 1 ' 2 ) < C^fi^Hy)^^ 

for all L and f,g 6 £ 2 (Al) with disjoint support. Moreover, all the quantities 

(6.3) E(\(q x q y )\), E(\(p xPy )\), E(| (q x q y ) p\ 1/2 ), E(\(p xPy )^ 2 ) 

satisfy bounds of the form Ce~^ x ~ y \ uniformly in L and x,y E A^. 

Note that, as indicated in Corollaries 14.51 and 15.61 the operators q x p y and p x q y have trivial 
ground and thermal state correlations, which we don't include here. 

Proof. These results follow by combining Corollaries I4.2[ 14. 5[ 15.31 and 15. 6[ whose assumptions 
hold in the Anderson regime due to Propositions IA.3T c) and IA.4f b). 

More precisely, the bounds on ground state correlations follow from Proposition IA.3( c). as 
the functions y?(s) = s 1//2 and (p(s) = s~ 1 ^ 2 have the required analyticity properties. Also, as 
y?(s) = coth(/3s 1//2 ) satisfies \(f(s) \ < C/((3s 1 ^ 2 ) for s near and has an analytic extension to the 
right half plane, the assumptions of Proposition IA.4( b) are satisfied with the choice tpi(s) = 
s^ 1 / 2 and ^2(3) = coth^s 1 / 2 ), leading to the decay bounds for thermal state correlations. 

□ 

As reviewed in Appendix |A], the Anderson model is localized at all energies, independent 
of the dimension, if the disorder is sufficiently large. This can be expressed by assuming that 
the density p of the distribution of the k x has sufficiently small L°°-norm ||p||oo ; reflecting the 
fact that the values of the random variables k x must be widely spread. A frequent way of 
referring to this in applications is by replacing the k x with Xk x , assuming that the distribution 
of k x has a density p as in A. 3 and that the disorder parameter A > is sufficiently large (the 
random variables Xk x have density p(-/A)/A which have small L°°-norm for large A). 

Applying Propositions IA.3I and IA.4I to the large disorder Anderson model leads to expo- 
nential decay of (—1/2, 1) as well as (—1, 1/2) eigenfunction correlators. This makes all our 
results on dynamics of disordered oscillator systems applicable, which we summarize in the 
following theorem. 

Theorem 6.2. Assume A. 3 with sufficiently large disorder, i.e. that \\p\\oo is sufficiently small. 
Then, in addition to the results of Theorem \6.1\ we have the following dynamical localization 
properties of Lieb-Robinson commutators: 
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There are constants C < oo and \J > such that 

(6.4) E ( sup || [r t (W(f)), W{g)\ \\] < C £ |/(x) | |y(y) |e^'l^l 

/or a// L and f, g : Al — > C, as well as 

(6.5) e( sup \\[T t (c x ),d y }\\ ) < C'e~^- y \ 



for all L and x,y G A i; where c x G {?x,Pa;} owe? d x G 

Moreover, the bounds on dynamic ground state correlations (4^4), j4-21 ) and jjj4-22 ) hold 



with r = 1. The bounds on dynamic thermal state correlations $5.10\) and $5.16}) hold with 



r = 1/2 and $5.17\) holds with r = 1 



Proof. In the large disorder regime we can combine Proposition lA.l( b) with Proposition lA.3( b) . 
choosing Eq = oo and a = —1/2, to show that h has exponentially decaying (—1/2,1)- 
eigenfunction correlators. Note here that the condition a > — s/{2 — s) of Proposition IA.3( b) 
is satisfied if s > 2/3 is chosen (which is possible as Proposition IA. 1( b) allows for s arbitrarily 
close to 1). 

That in the large disorder regime hi has exponentially decaying (— 1, l/2)-eigenfunction 
correlators follows from Proposition IA.4I 

The theorems in Sections [3l H] and [5] now yield all the claimed bounds. □ 

Another case in which localization can be proven at all energies is the one-dimensional 
regime, where no assumption on the disorder strength is required. The mechanisms behind 
one-dimensional localization are much more subtle than how localization follows (essentially 
by brute force) in the multi-dimensional large disorder regime. This leads to certain differ- 
ences in the mathematical characteristics of localization as summarized in Appendix [Al As a 
consequence we observe some slight changes in the consequences for disordered oscillator sys- 
tems, essentially in the form of "regularizing exponents". For simplicity, we state them only 
for our results on Lieb-Robinson commutators in the next result (which should be compared 
with (16.41) and (16.51) above) . Similar changes appear in some of the bounds on dynamical 
ground state and thermal state correlations. 

Theorem 6.3. Assume A. 3 and d = 1. Then there exist constants C < oo and fi' > such 
that 



(6.6) E(su P ||[t^(/)),^)]||) <C'J2\f(x)\ 1/2 \9(y)\ 1/2 e-»' lx 



for all L and f, g : Al — > C, as well as 

(6.7) E fsupmax {Hfnfe),^]!! 1 / 2 , || [^(p, ) , g y ] || , || [ n ( g;r ), p y ] || , || [n(p^), P y ] || }^) < C'e^'^ 
\teR J 

for all L and x, y G A^. 

Proof. The changes and slightly weaker results as compared to Theorem 16.21 are due to the 
fact that we now have to use part (c) of Proposition IA.1I which establishes localization for 
s-moments only for sufficiently small s > 0. This rules out an application of Proposi- 
tion |A]3](b) with a = —1/2. Instead we use Proposition IA.4( a). which gives exponential 
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decay of (— 1/2, l/2)-eigenfunction correlators. This yields (16. 6 p and the claimed bound on 
the first term on the left hand side of (16 .7p . 

For the other terms we can use Theorem l3.8( b). where exponential decay of (0, l)-correlators 
follows from Proposition IA.3( a) (which applies with arbitrarily small s). Without having 
stated it, we note that if either / and/or g are purely imaginary, then one could replace (16.61) 
by (16.41) . see Remark [3.41 above. 

□ 

We conclude this section with some remarks about possible applications of our results in 
Sections [3l H] and [5] to types of disorder other than the Anderson case discussed above, i.e. for 
other choices of the parameters m x , k x and X Xjy of the oscillator system. 

(i) Many (likely all) of our results will generalize to the case where the masses m x and 
couplings \ X)V are deterministic but non-constant, as long as they satisfy uniform bounds as 
in assumption A.l, i.e. < m min < m x < m max , < X XiV < A max , while disorder still enters in 
the form of i.i.d. random variables k x with assumptions as above. In fact, one can also relax 
the condition that the k x are identically distributed and merely needs that they are distributed 
with respect to densities which satisfy uniform bounds on their supports and L°°-norms. While 
we don't have explicit references for all the background results on Anderson localization in 
Appendix \A\ in this more general setting, it is well known and frequently discussed in the 
references that the existing tools allow such generalizations. 

(ii) Causing more technical effort, but still within the range of the existing tools to prove 
Anderson localization (while not considered in most of the cited references), is the case where 
the random variables k x have unbounded support, i.e. can be arbitrarily large. As a result, 
the operators hi do not satisfy norm bounds uniform in the disorder. In particular, (a, r)- 
eigenfunction correlators as defined in (I3.16P would now also have to be considered as singular 
for a = 1/2, just as the a = — 1/2 eigenfunction correlators are in the non-gapped case where 
the k x can be arbitrarily close to 0. 

(iii) More interesting would be cases where the k x are deterministic, for example a non-zero 
constant, while disorder enters in the form of random masses m x and/or couplings \ x , y . For 
suitable regimes, tools such as the fractional moments method used in Appendix |A] will likely 
yield localization results for the underlying single particle Hamiltonians hi with these types 
of random parameters. In fact, for the case of random X x ^ y some results on this are provided 
in Section 5 of [3]. 

But a major difference to the Anderson case is that these types of single particle Hamil- 
tonians hi will have an extended ground state. In fact, if k x = const, then hi as defined 
by (12. 6 p and (12.151) will have the ground state ip x = Cm x . As a result, one can not expect 
that hi is localized near the bottom of its spectrum. It is known in the theory of one-particle 
random Hamiltonians that extended states of this type, even if they only exist for a single 
energy level, may lead to non-trivial transport |24j . 

It is beyond the scope of our work to decide if this may lead to transport in oscillator systems, 
e.g. in the form of lower bounds on commutators as appearing in the Lieb-Robinson bounds. 
But we mention that the dispersive estimates proven in [5] for oscillator systems with constant 
parameters can be interpreted as a transport property inside the light cone established by the 
(non-zero velocity) Lieb-Robinson bound for such systems. We also mention that oscillator 
systems with random masses have seen considerable interest in the literature for many years 
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[38| |30| 112] . see [61 18] for some recent references. These works are concerned with classical 
oscillator systems and other types of phenomena than what we have considered here. 



In this section, we discuss possible extensions of the previous results to the infinite volume 
setting. For certain systems, the infinite volume single site Hamiltonian will have a positive 
lower bound uniform in the randomness. Under this assumption, the conditional results 
found in Sections [3j HI and |5] can be proven directly in the infinite volume. To see that this 
is true, one need only introduce the corresponding quantities of interest which, in this special 
case, are obviously well-defined, see e.g. the discussion below. A more serious challenge is 
to prove results that apply to gapless random oscillator models where the notion of singular 
eigenfunction correlators is not a-priori clear. With this as a goal, we focus solely on systems 
for which the single site Hamiltonian reduces to the Anderson model. In fact, throughout this 
section we will assume A. 3 as in Section [6] and therefore, the results found in Appendix IB] are 
applicable. 

We begin by introducing the CCR algebra on £ 2 (Z d ) and then define the time evolution 
corresponding to these harmonic models as a quasi-free dynamics on this Weyl algebra. A more 
thorough introduction to this general framework can be found in [TO]. To avoid additional 
complications, we focus on bounds for Weyl operators and only comment on how to treat 
positions and momenta at the end. 

For the Anderson model, it is proven in Appendix |Bl see specifically Lemma IB.3t that 
almost surely each 5 X is in the domain of the single site Hamiltonian. As a result, the notion 
of singular eigenfunction correlators (in this gapless case) is well defined and, e.g. in the large 
disorder regime, one can show that they decay appropriately. This is the main content of 
Appendix [B] Consequently, a zero velocity Lieb-Robinson bound, see Theorem 17.14 follows as 
in Section [3J 

The remainder of the section is used to introduce infinite volume ground and thermal states. 
These states are also defined in terms of singular objects which, in the specific gapless case 
where the single site Hamiltonian is the Anderson model, we can argue they are almost surely 
well-defined. Analogues of the results in Sections H] and |5] then follow. We focus our attention 
on thermal states as the corresponding results for ground states use similar arguments. 

We start by introducing the single-particle Hamiltonian h on l 2 {7L d \ Assuming A. 3, we 
extend the operator ho, whose matrix elements are still given by ( 12. 5p . to all of £ 2 {7L d \ The 
single particle Hamiltonian h is then defined as in ( I2.15P ; now on £ 2 {!/■). By assumption, h is 
self-adjoint and bounded; it still satisfies (I2.16P uniformly in the disorder. To avoid additional 
prefactors, we will also assume that the constant mass is chosen so that \i — I. 



For the remainder of this section, it will be convenient to identify £ 2 (Z d ) = £ 2 {7L d ; C) with 
£ 2 (Z d ; R) © £ 2 (Z d ; R) and regard 



With this identification, we can more easily express dependencies on real and imaginary parts 
in terms of matrix multiplication, see below. 

We now recall the definition of the Weyl algebra over £ 2 (/L d ). To each / e £ 2 (Z d ), associate 
a non-zero Weyl operator W(f) satisfying 



7. Results in Infinite Volume 



(7.1) 




as 




(7.2) 
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and 

(7.3) W(f)W(g) = e-^' 2 W{f + g) for all /, g G £ 2 (Z d ) , 
where 

(7.4) a(f, g) = lm[(f, g)] = (J/, g), J = ( J ~ J ) , 

and (-, •), resp. (■, ■), denotes the inner product on the complex (resp. real) £ 2 -space. It is well- 
known, see e.g. Theorem 5.2.8 of [TU], that there is a unique, up to ^-isomorphism, C*-algebra 
generated by the Weyl operators with the property that W(0) = I, W(f) is unitary for all 
/ G £ 2 (Z d ), and \\W(f) - I\\ = 2 for all / G £ 2 {Z d ) \ {0}. This algebra is called the CCR 
algebra, or the Weyl algebra, over £ 2 (Z d ) and we will denote it by W. 

A quasi-free dynamics on W, see e.g. Theorem 5.2.8 (4) [TU], is a one-parameter group of 
*-automorphisms r t of the form 

(7.5) r t (W(f)) = W(T t f) for all / G £ 2 {Z d ) 

where T t : l 2 (7L d ) — \ £ 2 (Z d ) is a group of real-linear, symplectic transformations, i.e. 

(7.6) (JTJ, T t g) = (/, g) for all f,g G f(Z d ) and t G R. 

Motivated by the finite-volume expressions (I3.14p and ( 13.15)) in Section [31 we define this 
infinite volume harmonic dynamics on W, as in (17. 5p above, in terms of the mapping T t given 
by 

cos(2t/i 1 / 2 ) -h 1 / 2 sin(2t/i 1 / 2 ) 
h~ 1 ' 2 sin(2t/i 1 /2) cos(2t/i 1 /2) 

where h is the infinite volume single-particle Hamiltonian previously discussed, and we recall 
that we have chosen constants so that \i = I. Since all entries in (17. 7p are bounded functions 
of h, T t defines a bounded, real-linear transformation on £ 2 (Z d ) for each t G K. by the usual 
functional calculus for self-adjoint operators. Moreover, it is easy to see that T t is a group 
and, one readily checks that 

(7.8) TjjT t = J, 

i.e. for each t G K, T t is symplectic. As is indicated above, this is sufficient to define the 
harmonic evolution as a quasi-free dynamics on W. 

We can now prove the analogue of Theorem 13.31 in this setting. Recall that in the notation 
of A. 3, p is the density of the random variables and that large disorder is expressed in terms 
of smallness of 



e ~H\x-y\ 



(7.7) T t 



Theorem 7.1. Let ||p||oo be sufficiently small. Then there exist C < oo and fi > such that 
(7-9) E(sup\\[r t (W(f)),W(g))\\] < C^2\f{x)\\g(y)\ 

for all f,ge £ 2 (Z d ). 

Proof. As is discussed in the beginning of Appendix [Bj an infinite volume version of Propo- 
sition |Aj] is well-known to hold. Thus for any d > 1 and \\p\\oo sufficiently small, Propo- 
sition IA.l( b) holds and so the assumption of Proposition IB.4( a) is satisfied on an interval 
containing the (uniformly bounded) spectrum of h. In this case, the proof of (17. 9p immedi- 
ately follows using the arguments of Theorem 13.31 with ( IB. 12[) as the new input. □ 
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We now turn our attention to states and correlation bounds. To begin, we review how 
these infinite volume states are rigorously introduced. We first sketch the details for thermal 
states. The corresponding results for ground states immediately follow. Again, more general 
information, including these facts, can be found e.g. in [TO] . 

Fix (3 > and consider the functional on W defined by setting 

(7.10) M W (f)) = e-3**C/./) for all / G l 2 (Z d ) 

where sp is the quadratic form given by 

(711) , (f n _f \\WhfV for all /eP 

[< ' W>/;-| +oc otherwise 

with 

/ h-V* 00^(^)1/2 

and form domain 

(7.13) ^ = {/ € f(Z d ) : Re[/] G D {h' 1 ' 2 )} , 

where D{h~ 1 / 2 ) is the domain of /i" 1 / 2 , see e.g. Lemma IB. 21 and the preceding discussion. 
Note that if / G then /i" 1 / 4 coth(/3/i 1 / 2 ) 1 / 2 Re[/] is well-defined since h 1 ' 4 coth(/3/ i 1 / 2 ) 1 / 2 
is bounded and all these functions of h commute. Since h is random, so too are ujp and its 
domain Dp, however, given Lemma IB. 31 we know that almost surely J 7 , the (deterministic) 
set of all functions in £ 2 (7j d ) with finite support, is in Dp. 

To see that the functional up actually defines a state on W, i.e. a positive linear functional, 
we need only check that 

(7-14) a(f,g) 2 <Asp(fJ)sp(g,g) 
see e.g. [13]. This follows by noting that the matrix 

h- 1 '* 

h l ' A 



(7.15) M 



is symplectic, that ||M/|| < H-M^/H for all / £Dp, and therefore 

(7.16) a(f,g) 2 = a(Mf,Mg) 2 < \\M f\\ 2 \\Mg\\ 2 < Asp{f,f)sp{g,g) . 
We next note that up is time invariant, i.e., the fact that 

(7.17) up(r t (W(f))) = up(W(f)) forall/G£ 2 (Z d ). 

First, observe that if / G Dp, then T t f G Dp for all t G K. This follows by inspection of 
the mapping T t , see (ITTTjl . Since T t is invertible, this verifies (I7T7I) for all / G £ 2 (Z d ) \ Dp. 



Otherwise, a short calculation shows that 



cos(2^V2) _ sin (2t/ i i/2> 



1/2^ 



(7.18) M s r, = fl«M s where ft = (^ 2(A1/2) cos(2i;i 

which is clearly orthogonal for all i e R. (17.171) is proven. 

We can now state our results on dynamic correlations for the infinite volume thermal states. 
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Theorem 7.2. Let \\p\\oo be sufficiently small. Then there exist C < oo and \i > such that 

E (sup \u p (r t (W(f))W(g)) - u p (r t (W(f))) u p (W(g))\ 

(7.19) Ktm 

<cY,\f^)\ ll2 ^v)\ ll2 ^ x ^ 

for all f,g EJ 7 . 

Proof. As in the previous theorem, it is clear that for ||p||oo sufficiently small, the assumptions 
of Proposition IB. 5( a) hold on an interval containing the spectrum of h. In this case, one sees 
that the argument for the proof of Theorem 15.11 (with r = 1/2) carries through, using the 
results of the above discussion and (IB.23[) as input, to show (I7.19p . □ 

In addition, we also have the follow result on static correlations, which does not require the 
large-disorder assumption but only depends on localization near the bottom of the spectrum 
of the single-particle Hamiltonian. 

Theorem 7.3. There exist C < oo and fi > such that 

(7.20) E (\ Up (W(f)W(g)) - Ufi {W{f))u p (W(g))\) < ^{v)^-^ 
holds for all f,gEJ- with disjoint support. 

Proof. One easily checks that an analogue of Corollary 15 .31 holds. Now, for these static bounds, 
we use Proposition IB. 5( b) with (fi(h) = <^>2(h) = h -1 ! 2 . It is now clear that (I7.20p follows 
from the analogue of (15. lip using (1B.24|) . □ 

Similar results, i.e. analogues of Theorem I7.2l and Theorem 17. 3| hold for ground states. Note 
that ground states can be defined as the /3 — >■ oo limit of thermal states, or, more concretely, 
by setting coth(/3/i 1//2 ) = 1 in (I7.12p . For dynamic and static correlation estimates in this less 
singular state, we use Proposition IB.4( a) and (b) respectively. We leave the details to the 
reader. 

The ground state and thermal states of the infinite systems considered here are sufficiently 
regular (see [10j pp 37-38]), so that one can define position and momentum operators in the 
GNS representation almost surely. As a consequence, one can verify that, on a dense subset 
of the GNS Hilbert space, formulae identical to those in Lemma [3.71 hold. On this subset, an 
analogue of Theorem I3.8l can be proven, again under the large-disorder assumption. Analogues 
of the results for correlations of position and momentum operators in the ground state and 
thermal state similarly hold in the GNS representations. 



8. Conclusion 

We conclude with a brief summary of what we consider to be the most important novel 
contributions of our work. 

We prove dynamical localization in the form of a zero-velocity Lieb-Robinson bound for 
harmonic oscillator systems. Disorder provides the crucial mechanism for this, as it is necessary 
that the associated one-particle Hamiltonian is localized at all energies. 

We prove exponential decay for ground state and thermal state correlations of disordered 
oscillator systems. For this it is not required that the oscillator system is gapped (as in related 



30 



B. NACHTERGAELE, R. SIMS, AND G. STOLZ 



previous results for deterministic systems). Instead we exploit a "mobility gap" expressed in 
terms of localization of the one-particle Hamiltonian near the bottom of its spectrum. 

In cases where the mobility gap extends to all energies, as for the dynamical localization 
bounds, we can go beyond static correlations and get exponential clustering for dynamic 
ground and thermal state correlations. 

The extensions of our results to oscillator systems over the full Euclidean lattice provide 
localization results for oscillator systems in infinite volume. 

In order to cover non-gapped as well as infinite volume oscillator systems, we need new 
results on Anderson localization. Specifically, we prove exponential decay of some types of 
singular eigenfunction correlators for the Anderson model required for our applications. 

Appendix A. Some Results on Anderson Localization 
We consider the finite- volume Anderson model on £ 2 (Al), Al = [-L, L] d D Z d , 
(A.l) h L = h 0lL + V u . 

Here fro,£ is the graph Laplacian on A^, given in terms of its quadratic form by 
(A-2) (f,h , L g)= (f(y)-f(x))(g(y)-g(x)), 

{x,y}CA L 
\x— 5/1=1 

and 

(A.3) (VJ)(x)=cu x f(x) 

for an array (uj x ) x€ z d of i.i.d. random variables. We assume that their common distribution 
is absolutely continuous with bounded and compactly supported density p. We will assume 
min(supp p) = 0. Thus for any given x the event u x = has probability zero. This implies 
that almost surely hi > (i.e. does not have eigenvalue 0) and thus is invertible. Nevertheless, 
the bottom of the almost sure spectrum of the infinite volume Anderson model h + (h 
the graph Laplacian on Z d ) is 0, see [32] for an elementary proof. 

This latter choice of normalization as well as the choice of the graph Laplacian ( 1A.2|) as 



background operator are motivated by our applications in earlier sections. Other background 
operators could be used for all the localization results stated below without affecting the 
proofs. In particular, as done in most of the references given below, one could work with 
the restriction Xa l ^oXa l of the infinite volume Laplacian to A^, which differs from /i ,l by a 
boundary condition (/io,L is sometimes called the discrete Neumann Laplacian) and an energy 
shift (if h is chosen as the next-neighbor- hopping or adjacency operator). We could also use 
periodic boundary conditions in (IA.2j) . i.e. replace by a <i-dimensional discrete torus. In 



the latter case distances of lattice sites in flA.2j) as well as in the results below have to be 
interpreted as distances on the torus. 

The main reason for requiring that the distribution of the u x has bounded and compactly 
supported density p is that this allows us to refer to the rather strong localization properties 
of the Anderson model which have been proven by the fractional moment method (FMM), 
see [21 E [13 H] and, for a recent survey, [32]. Some of the results proven by this method hold 
under the weaker assumption of Holder continuous distributions, e.g. [3], and the assumption 
of compact support of p could also be relaxed. But the method can not be extended to more 
singular distributions and other methods do not yield equally strong results. 
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The main technical characteristic of the FMM is that it initially aims at proving fractional 
moment localization, from which other properties such as spectral and dynamical localization 
are derived as consequences. 

For given < s < 1 and an interval / C t, we say that the finite-volume Anderson model 
hi has localized s-moments in / if there are constants C < oo and p > such that 

(A.4) E(\G L (x,y;E + ie)\ s ) < C'e'^'^ 

for all positive integers L, all x, y G Al, E E I and £6R. Here Gl(x, y\ z) := (5 X , (Jil — z)~ 1 5 y ) 
is the Green function for Hl and E(-) refers to the disorder average. 

The following summarizes the regimes in which fractional moment localization for the finite 
volume Anderson model has been shown. 

Proposition A.l. The d- dimensional finite volume Anderson model hi has localized s- 
moments in I under each of the following assumptions: 

(a) d > 1 and s G (0, 1) arbitrary, and I = [0, E ] for some E Q = E (s, d) > 0, 

(b) d > 1 and s G (0, 1) arbitrary, WpW^ < p max for p max = p max (s, d) > sufficiently small, 
and I = [0, oo) , 

(c) d = 1, s > sufficiently small, and I = [0, oo). 

Historically, the first localization result obtained by using the FMM was for the large dis- 
order regime (b), see [3] and [TJJ]. In particular, the statement in (b) covers the case of an 
Anderson model /io,L + AV^, where is given via a fixed density p and the coupling parameter 
A is sufficiently large. Band edge fractional moment localization (a) was obtained in [1]. That 
(c) follows from the well known properties of the Lyapunov exponents of the one- dimensional 
Anderson model was first observed in [32], see also [20]. Much of this, including the proofs, is 
reviewed in [12]. 

One of the virtues of working in finite volume is that e = is allowed in (1A.4|) . The proofs 



show that any fixed E G M is an eigenvalue of with probability zero and thus, as the 
spectrum of is discrete, (hi — E)^ 1 exists almost surely and the left hand side of (IA.4f) 
makes sense for e = 0. 

Among the consequences of (IA.4[) . the one of most interest to us here is the following: We 
say that the Anderson model hi has localized eigenfunction correlators in an interval let, 
if there exist p > and C < oo such that 

(A.5) E ( sup | (5 X , u(h L ) X i(h L )5 y ) | ) < Ce~^ 

\H<i J 

for all positive integers L and x, y G A^. The supremum is taken over all functions u : M. — > C 
which satisfy the pointwise bound |tt| < 1 on J. The operators u{hi) and xK^l) are defined via 
the functional calculus for self-adjoint operators, in particular Xi^l) is the spectral projection 
for hi onto the interval /. Of course, as we work in finite volume, the functional calculus 
comes down to simple finite eigenfunction expansions. For this reason we also don't need any 
measurability assumptions for the functions u. 

For more explanation of the reason for describing the left hand side of fl A.5[) as an eigen- 
function correlator see [H IU H2]. The choice u t (E) = e~ ltE , t G M, shows that (1A.5I) implies 
the dynamical localization bound 

(A.6) E (sn V \(S x ,e- uhL xi(h L )S y )\) < Ce'^. 
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While this explains why OA.Sp is sometimes referred to as dynamical localization, we will need 
(IA.5P also for other types of functions u. 

The crucial bound which allows to deduce localization of eigenfunction correlators from 
fractional moment localization is the content of the next result. 

Proposition A. 2. For every < s < 1, every dimension d, and every bounded and compactly 
supported probability density p, there exists C" = C"(s,d,p) < oo such that 

(A.7) E (supJ(S x ,u(h L ) X i(h L )Sy)\\ < C" ^jE(\G L (x, y] E)\ s )dE 

for every interval Id. 

This relation originates from the works [U 0], a detailed proof can also be found in [12]. By 
combining Propositions IA.l1 and IA.21 one gets the first part of the following result. 



Proposition A. 3. Let ht be the finite-volume Anderson model defined by liA.l]) . liA.ty) and 
liA.3\) with assumptions as above. 

(a) If hi has localized s-moments in I in the sense of ( A. 4% then 



(A.8) E sup \(S x ,u(h L ) X j(h L )5 y )\ < C\ j| 1/(2-^-/^1/(2-*) 

for every subinterval J of I, every positive integer L and all x, y E A^. Here C = C" / (C / ) 1 /( 2 ~ S ) 
with C from ( L4-4 > and C" from \A.T\} . 



(b) If hi is s-localized in [0, Eq] in the sense of ( A. 4 ) and a > —s/{2 — s), then there exists 
C\ < 00 such that 

(A.9) E ( sup \(S x ,hlu(h L )xioM h L)S v )\) < C ie -^-y^ 2 -^ 

\H<i J 

for all L and x,y G Al. 

(c) Suppose that (p : (0, 00) — > C satisfies \<p(t)\ < Ct a for t near with some C < 00 and 
a > —1, and that p> has an analytic extension to a semi-strip {z : Rez > 0, \Imz\ < 77} for 
some rj > 0. Then there exist Ci < 00 and P2 > such that 

(A.10) E(\{5 x ,ip(h L )5 y )\) <C 2 e~^ lx ~ yl 

for all L and j,?/6Aj,. 

The statements in (b) and (c) involve the terms hi and <p{hi) which are allowed to be 
singular at energy E — 0. Also, the left hand side of (c) does not require to include the 
projection X[o,e ](^) on ^° localized energy regime as in (b). The price for this is that one 
can not include a supremum over a class of functions of hi as in (a) and (b) and, in particular, 
say nothing about dynamics. 

Below we prove (b) using (a) and a Riemann sum argument. The proof of (c) is done by 
contour integration, similar to methods previously used in [2]. 

Proof of Propositon \A . 3\ It remains to prove (b) and (c). 

If a > 0, then fl A. 9j) follows easily from (IA.8[) with J = [0, E ], as in this case x a < E$ on 
[0,E ]. 
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For — s/(2 — s) < a < 0, decompose the interval (0,E ] into Z„ := (E /(n + l),_E /ra], 
n = 1, 2, . . .. As is almost surely not an eigenvalue of h^, we have with probability one that 

oo 

(A- 11 ) X[o,E a ](h L ) = X(o,E ](h L ) = ^2xi n ( h L)- 

n=l 

Convergence is trivial here as, due to the discreteness of the spectrum of h^, the sum is finite. 
If \u\ < 1, then < (E /(n + 1)) Q on I n . Thus it follows from (a) that 

(A.12) E (sup \{S x ,hlu(h L ) X[ o,E o] (hL)S y )\) < £)<7 ( -^Y |j B |V<^) e -^l/P--). 
\H<i / ^ V^ + V 



As (-^j-)" |4| 1/(2 s) ~ (-)" 2 ~ s and a + > 1, the series converges. We get flA.9j) with 



/(2- S )_ 

To prove (c), we first choose s < 1 such that a > — s/(2 — s) and then i?o > as in 
Proposition IA.l( a) . Decompose 

(A.13) E(\(5 x ,<f(h L )S y )\) < E(\(S x MhL)X[oM h L)S v )\) 

+ E(\(5 x ,(f(h L )x(E ,oc)(h L )5 y }\) . 

For the first term we get the required bound from (b), using the bound \<p(t)\ < Ct a . For 
the second term we argue as follows: 

As p is compactly supported there exists M < oo such that < M uniformly in L and 
the disorder. Assume that E is not an eigenvalue of (which holds almost surely). Then, 
with E(-) denoting the spectral resolution of 



(A.14) {5 x ,ip(h L )x(E ,oo)(h L )5y) = / <p(t) d(5 x ,E(t)5 y ) 

>(E ,M] 



I 7T-- f 2&dzd(6 x ,E(t)6y) 

J(E ,M] 27TZ J T Z-t 



by the Cauchy integral formula, where we have used that if has an analytic extension to a 
semi-strip and thus we can choose for T the rectangular contour with vertices E — ir]/2, 
M + 1 - ir)/2, M + 1 + ir)/2 and E + irj/2. 

We work in finite volume, so that the spectral integral in (IA.14I) is a finite sum, allowing 
for a trivial exchange of integration order. One gets 



(A.15) \(S x ,<f(h L )x( Eo ,oo)(hL)Sy)\ 



2m 



r 



<p(z){6 x , (h L - z) 1 X(E ,M]( h L)Sy) dz 



< C \(6 x ,(h L -z) 1 X(E ,M]( h L) s y)\\ dz \ 



with C = w- max.{\tp{z)\ : z G T}. 

Decompose T = Fi U T2, where Ti is the line segment from Eq + if]/2 to Eq — ir]/2 and 
T2 the remaining part of the contour. Choose < s' < 1/2 and use the trivial bound 
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\\{h L — (Eq + < 1/a to estimate 

(A.16) 

E / \{5 x ,{h L - zY l X{EoM] {h L )5 y )\\dz\ 

rri/2 
V /2 

n/2 



< / a~^E(\(6 xl (h L -(E + ta))-\ (EoM] (h L )5 v }\ s ') da 
J -n/2 v J 

< J V/ \~^E f^|G L (x,^Eo + m)| s '|(^,x (So ,A/](/iL)^)| s ' J da 



-n/2 

< / V{(E|G L (x, V ;^ + »a)| 2 'y/ 2 (E|((y w ,X(Eb,^(^K>r s ? /2 } da. 

By Proposition ED E(|G L (x, E + ia)| 2s ') < Cxe^ 1 ^!, while 
(A.17) E(|(^ )X ( fi0 ,M](^)^}| 2s ') < E(\(5 v ,X(E ,M](hL)S y )\) 

= E(\(5 v ,(I-X [0 ,E 0] ( h L))Sy)\ 

< S vy + C 2 e-^ v - y \ < (C 2 + l) e -^l^l 
by (a). We can thus bound the right hand side of (IA.16[) by 

(A.18) C a -(i-«)y^ e -Mil*-«l/2 e -w|«-tfl/2 < c^-Mls-wl 

for any /i < | minj/i!, /i 2 }. A similar bound is found for the second term contributing to 
( IA.15j) . corresponding to r 2 , where one can use the well-known deterministic Combes-Thomas 
bound \Gl(x,v; z)\ < C\e~^ l][X ~ v \ which holds uniformly in z G T 2 . For a self-contained proof 
of the latter, which extends to the finite volume case considered here, see, e.g., Section 11.2 

of [25]. 

All of this can now be combined to get the desired decay bound for the second term in 
f lA~T3|) . 

□ 

Parts (b) and (c) of Proposition IA.3I exclude the case a = — 1. As our final result in this 
appendix we show how these stronger zero-energy singularities can be handled by including a 
fractional moment (the 1/2 moment) in the expectation of the eigenfunction correlator. 

Proposition A. 4. Fix assumptions as above. 

(a) If h E has s-localized moments in [0, Eq] and a > — 1 — ^f-, then there exist C3 < 00 and 
/X3 > such that 

(A.19) E ( sup h a L u(h L ) X[ o,E ](hL)Sy)\ 1/2 j < 

for all L and x, y G A^. 

(b) Suppose that tpi and <p 2 both have the properties of the function tp in Proposition \A.3[ 
Then there exist C4 < 00 and \i± > such that 

(A.20) E (|(4,^ 1 (/ iL )^ 2 (/ iL ) ( 5A| 1 /2) < CAe -^-v\ 

for all L and x, y G A^. 
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Proof. Choose ot\ > — 1 and «2 > — s /(2 — s) such that a = a.\ + «2- Then 
(A.21) (5,,/i^(/i L )x[o^ ](^)^} = (h^8 x ,hfu{h L ) X \o,E ){h L )8 y ) 

= ^2(hT5 x ,5 w }(5 w ,hl 2 u(h L )x[o,Eo](h L )5 y ). 

w 

Using Cauchy-Schwarz on the expectation we find from this that 
(A.22) E ( sup \(5 x ,h a L u(h L ) X[ o,E ]( h L)Sy) ll/2 

\\u\<l 

J2m(^h a L ^ w }\)) 1/2 ( E(sup \(5 w ,h a L *u(h L ) X[ o,E o] (h L )5 y )\) 1 

l«l<i 

Part (c) of Proposition IA.3I applies to the first term of the sum, while part (b) applies to the 
second term, giving exponentially decaying bounds in \x — w\ and \w — y\, respectively, which 
carry over with a suitably reduced rate to the convolution-type sum. This proves (a). 
For (b) one argues similarly, bounding 

(A.23) E(\(5 x ,vi(h L )MhL)5 y )\ l/2 ) 

< ^(E(|(^, V 9 1 (/ iL )^)|)) 1/2 (E(|(^,^(/ ii )^)|)) 1/2 , 

w 

and using Proposition IA.3( c) on both terms to conclude as above. 

□ 



Appendix B. The Anderson model in infinite volume 

Our goal here is to extend the results of Appendix |A] to the Anderson model in infinite 
volume. On £ 2 (Z d ) consider 

(B.l) h = ho + V u . 

The graph Laplacian ho on Z rf is given by 

(B.2) (f,hog)= E (f(v)-f(x))(9(v)-9(x)), 

{i,j)CZ <! 

\x— y|=i 

and the random potential is defined by fl A. 3j) with i.i.d. random variables (u x ) xe %d. As before 
we assume that they are distributed according to a bounded and compactly supported density 
p with min(suppp) = uj min > 0. By well known facts (e.g. [T7]) the spectrum of h is almost 
surely given by [0, 4d] + supp p, in particular, miner (/i) = w m i n , and E = u m i n is almost surely 
not an eigenvalue of h. 

We start by noting that Proposition IA. II extends to infinite volume. By this we mean that 
the infinite volume Green function G(x, y; z) := (5 X , (h — z)~ l 5 y ) satisfies a fractional moment 
bound of the type (1A.4|) for all x, y G Z rf and uniformly in e ^ in the same regimes as 
described in Proposition IA.1I This can be found in the same references as given for the finite 
volume case. Note that in infinite volume one can not allow e = as real energies will almost 
surely be contained in the spectrum of h (in situations other than trivial ones). 

Another known fact is that Proposition IA.3( a) extends to eigenfunction correlators for the 
infinite volume Anderson model. 
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Proposition B.l. Assume that hi has localized s-moments in an interval I for some s 6 
(0,1). Then 

(B.3) E( sup \(5 x ,u(h)xAh)5 y )\) < C\ j|i/(2-») e -Ml*-wl/(2-.) 

V m:R->C Borcl, |u|<l / 

for every bounded open subinterval J of I and all x,y G Z d . Here C and \x are the constants 
from Proposition \ A. 3( a). 

That ( IB. 31) follows by taking the L — > oo limit in (1A.8I) is essentially a consequence of 
strong resolvent convergence, see [31 [32] for details. The restriction to Borel functions u and 
open intervals J is due to measure theoretic requirements of the functional calculus and, in 
particular, the use of Lusin's theorem in the proof. 

All the remaining results in this Appendix are most interesting for the case u m \ n = (for 
Wmin > our results are obvious or easier to prove). This is the most interesting case for our 
applications, as it corresponds to gapless oscillator systems. It is also the most difficult case, 
as in this case operators such as /i -1 / 2 , used to define singular eigenfunction correlators, are 
almost surely unbounded. We start by stating a lemma to clarify the meaning of h~ l l 2 in this 
case. 

For this we may consider any non-negative self-adjoint operator h in a separable Hilbert 
space T-L such that is not an eigenvalue of h. Thus h : D{h) — > R(h) is injective and 
h^ 1 : R(h) — > D(h) is a well defined (but generally not bounded) linear operator. One may 
also express h~ x in terms of the functional calculus of self-adjoint operators. For this let E(-) 
be the spectral resolution of h and let <p : R — > R be defined by 



(B.4) <p(t) 



1/t, if t > 0, 
0, if t < 0. 



Then a self-adjoint operator <p(h) in "% is defined by the functional calculus (e.g. Chapter 7 
of EH) as 



(B.5) D{<p{h)) = {f eV.: f ^d\\E(t)f\\ 2 < oo}, <p{h)f = / ~dE{t)f. 

</(0,oo) ^ J(0,oo) £ 

Lemma B.2. (a) h" 1 = (f(h). 

(b) (p(h) is a non-negative self-adjoint operator and its unique non-negative square root 
((^(/i)) 1 / 2 is given by ij)(h), where 

r 1 / 2 , ift>0, 
0, ift<0. 



(b.6) m 



(c) The unique non-negative square root h 1 / 2 of h is invertible and ip(h) = (h 1 ^ 2 ) 1 . 

Proof. We freely use background from the functional calculus of unbounded self-adjoint oper- 
ators, see e.g. Section 7.2 of [33] (note that the notation E(u) used there corresponds to u{h) 
in the case where E(-) is the spectral family associated with h). 

(a) The functional calculus allows to write h = id(h) with the identity function id(t) = t for 
all i e R. We have <^-id = X{o,<x>) an d thus (tp ■ id)(h) = X(p,oo){h) = h as h > and is not an 
eigenvalue of h. By Theorem 7.14(h) in [44J we thus get (p(h)h C / and D(ip(h)h) = D(h). The 
first fact implies <p(h) C far 1 , while the second gives D^h" 1 ) = R(h) C D(ip(h)). Combined 
this yields the claim. 
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(b) ip > implies <p(ti) > by Theorem 7.14(f) of [H]. It remains to show that ip(h) 2 = 
<p(h). To this end, we first observe that D(<p(h)) C D(tp(h)). This is the same as saying 
that f^^t' 2 d\\E(t)f\\ 2 < oo implies f^ Q ^t' 1 d\\E(t)f\\ 2 < oo, which follows from Cauchy- 

Schwarz. Now, as ip 2 = <p, it follows from Theorem 7.14(h) in [H] that ip{h) 2 C <p(h) with 
D(ip(h) 2 ) = D(i)(h)) n D(<p(h)) = D(<p(h)), giving the claim. 

(c) This follows from an application of Theorem 7.14(f) of [H] similar to the one in part 
(a), this time to the product of the functions ip and g, where g(t) = t 1 ^ 2 for t > and g(t) =0 
for t < 0. 

□ 

Based on these facts we will from now on write h^ 1 ^ 2 for ip(h), which means that we have 
defined h~ x l 2 = (/i -1 ) 1 / 2 . This notation is further justified by Lemma IB. 2( c). which says that 

h -l/2 = (^l/2)-l_ 

We now return to studying the infinite volume Anderson model h. 

Lemma B.3. Let h be the Anderson model on 7L d with assumptions as above. Also let x G Z d . 
Then 5 X E D{h~ 1 ^ 2 ) almost surely. 

Proof. By Lemma TB. 2( b) we have to show that C 0oo > t~ x d\\E{t)b~ x \\ 2 < oo almost surely. This 
will follow from 

(B.7) Elf/" r 1 d\\E(t)5 x \A ' ) <oo. 



y(o,oo) / 
By the infinite volume version of Proposition IA. 1( a) for s = 1/2 we have 

(B.8) E^M^-^rW 72 ) <C<oo 

uniformly in e > 0. 

Using that almost surely is not an eigenvalue of h, we get that, almost surely, 
(B.9) ReiS^ih-ieyX) = Re / ^-d\\E(t)5 x \\ 2 

J(O.oa) t — 1>£ 



d\\E(t)5 x "' 



(0,oo) 1 + E 

For all t > we have t/(t 2 + e 2 ) / 1/t as e \ 0. Thus the monotone convergence theorem 
shows that 



2 

a; II ) 



(B.10) / -* dll^kHV / ^II^H 

J(P,oo) 1 ' £ J{0,oo) 1 

as well as 

(bu) ^{(l^ B{t ^T) = ^{(I^^^T 

Combining (1B.8I) and (1B.9I) (and |Rez| < \z\) we see that the right hand side of (IB. Ill) is 
uniformly bounded in e > 0. This yields (IB.7p . 

□ 

We can now state and prove infinite volume versions of parts (b) and (c) of Proposition IA.3I 
While more general results hold, we only do this for a = —1/2, the relevant case for our 
applications. 
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Proposition B.4. (a) If Hl is s-localized in [0,-Eo] for some s G (2/3,1) , then there exists 
C\ < oo such that 

(B.12) E( sup \{5 x ^ l l 2 u(h) X[W (h)5 v )\ ) < c^-mI^I/M 

\u Borel, |u|<l J 

for all x, y G Z d . Here fi > is as in Proposition \B.l[ 

(b) Suppose that ip : (0, oo) — > C satisfies cp(t) < Ct~ x l 2 for t near 0, and that ip has an 
analytic extension to a semi strip {z : Re 2 > 0, |Im£| < r]} for some f] > 0. Then there exist 
C*2 < oo and fi2 > such that 

(B.13) E(\{6 x ,<p(h)6 y )\)<C 2 e-^ x -y\ 
for all x,y G Z d . 

Proof. As is almost surely not an eigenvalue of h and h~ 1 ^ 2 commutes with the bounded oper- 
ator u(h)x[o,E ](h), the vector h~ l l 2 u(h)x{o,Eo\{h)5 y is almost surely well-defined by Lemma [E~3l 
Thus the left hand side of (IB. 121) makes sense. The same can be said for the left hand side of 
(IB. 131) as we can write (p(h) = h~ x l 2 h}l 2 Lp{K) and, by assumption, h}l 2 Lp{K) is bounded. 

The proof of (a) now proceeds by essentially the same Riemann sum argument as the proof 
of Proposition IA.3( b). We use Proposition IB. II instead of Proposition IA.3( a). which requires 
to decompose [0, E Q ] into open intervals /„ := (E /(n + l),Eo/n). This does not impact the 
proof as the countable set {0}U {E /n : n — 1, 2, . . .} almost surely carries no spectral measure 
of h. The only other part of the proof which needs a bit more care than before is flA.lip . Here 
we use again that almost surely 5 y G D{h~ l l 2 ). Thus, almost surely, 

(B.14) h-VMQxvwiWv = J2 h ~ 1/2u ( h ^W 6 y 

in the sense of strong convergence. The rest of the argument is unchanged. 

To prove (b), we first note that we may assume that \(f(t)\ < Ct~ x l 2 not just near 0, but 
uniformly on all the spectra of h and hi, which is due to the uniform boundedness of these 
spectra. 

Assume that is not an eigenvalue of h and that 5 y G D^h^ 1 ^ 2 ), which, by Lemma IB. 31 
holds almost surely. By Lemma IB. 21 this means f^ ^t -1 d\\E(t)5 y \\ 2 < oo. It follows that 

J (0oo) I^I^Hi^cy 2 < oo, i.e. 5 y G D(<p(h)). 

For 5 > define continuous functions <p$ ■ M — > C such that tps{t) = (p(t + 5) for t > (this 
is possible as 5 > allows to avoid the singularity of ip at 0, while the exact choice of ips on 
(—00,0) is irrelevant). 

We claim that 

(B.15) <f{h)5 y = lim <p 1/n {h)8 y . 

n—^oo 

To see this, note that by the functional calculus 



(B.16) \\{<p(h) - <p 1/n (h))6 y \\ 2 = / \<p(t)-ip(t + l/n)\ 2 d\\E(t)5 y \\ 2 . 

</(0,oo) 

We have tp(t + 1/n) — > (p(t) for all t > by continuity as well as 

AC 2 

(B.17) \<p(t) - <p{t + l/n)\ 2 < 2(\ V (t)\ 2 + \<p(t + l/n)\ 2 ) < 



t 
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on a(h). The latter is integrable with respect to d\\E(t)5 y || 2 . Thus (IB.15P follows by dominated 
convergence. 

Also, for fixed n the function (py n is continuous and thus it follows from strong resolvent 
convergence that 

(B.18) ip 1/n (h)5 y = lim (pi/ n {h L )5 y . 

It remains to show the existence of C 2 < oo and /x 2 > such that 
(B.19) E (\(5 x ,<f 1/n (h L )5 y )\) < C 2 e~^ x - y \ 

for all i,t/e Z d , uniformly in L and n. For (IB.15P and flB. 18j) combine to give 
(B.20) \(6 x ,<p(h)6 y )\ = lim lim \(6 X , <p 1/n (h L )5 y )\ 

n— >oo L— >oo 

almost surely. Thus (1B.13P follows from (IB.19P and Fatou's lemma. 

The proof of (1B.19P is done by a straightforward inspection of the proof of Proposition lA.3( c). 
showing that one may replace (p by <pi/ n and gets bounds uniform in n. This comes down to 
two instants in the proof. In the part corresponding to the spectral projection on [0, Eq] use 
that 

(B.21) fi/n(h L )x[o,E ](h L ) = h L ^n{hL)X[o,Eo](h L ), 

where \ip n (t)\ = \t 1/2 <pi, n (t)\ < C. Thus one can use Proposition IA.3( b) uniformly in n. 
Second, when bounding E(|(<J X , </ ? i/n.(^i)X[s ,oo)(^L)^)|) all one needs to observe is that the 
constants 

(B.22) C n = ^-m a x{\ Vl/n {z)\:zeT} 

entering (IA.15I) are uniformly bounded in n. 

□ 

It is now straightforward to adapt Proposition IA.4I and its proof to infinite volume. All that 
needs to be done to ensure that the infinite volume quantities are well-defined is to distribute 
singular operator functions evenly to both sides of the inner products. 

Proposition B.5. (a) If hi has s-localized moments in [0, -Bo] f or some s £ (2/3,1), then 
there exist C 3 < oo and fi 3 > such that 

(B.23) E( sup K/i- 1/2 ^,/i- 1/2 «(/i)X[o,^](/i)^)r /2 J <C z e-^ x ~y\ 

\u Borel, |u|<l I 

for all x, y £ Z d . 



(b) Suppose that tpi and ip2 have the properties of ip in Proposition B.^ (b). Then there exist 
C4 < oo and n± > such that 

(B.24) E (\(Mh)5 X iMh)5 y )\ l/2 ) < 

for all x, y £ Z d . 
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